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FROM NEWTON TO BOLTZMANN: THE CASE OF SHORT-RANGE 

POTENTIALS 



Isabelle Gallagher, Laure Saint-Raymond, Benjamin Texier 



Abstract. — We fill in all details in the proof of Lanford's theorem. This provides a rigorous deriva- 
tion of the Boltzmann equation as the mesoscopic limit of systems of Newtonian particles interacting 
via a short-range potential, as the number of particles N goes to infinity and the characteristic length 
of interaction e simultaneously goes to 0, in the Boltzmann-Grad scaling Ne^ 1 = 1. The case of 
localized elastic interactions, i.e., hard spheres, is a corollary of the proof. The time of validity of the 
convergence is a fraction of the mean free time between two collisions, due to a limitation of the time 
on which one can prove the existence of the BBGKY and Boltzmann hierarchies. Our proof relies on 
the important contributions of King, Cercignani, Illner and Pulvirenti, and Cercignani, Gerasimcnko 
and Pctrina. 
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CHAPTER 1 



THE BOLTZMANN-GRAD LIMIT 



We study the qualitative behavior of systems of interacting particles of the form 

dx t dvi > 

(l.o.i) ~dt= Vl ' m-fir = - 2^ V&(xi - xj) , 

for 1 < i < N, where (xi, Vi) € R d x R d denote position and velocity of particle i with mass raj (which 
we shall assume equal to 1 to simplify) and the force exerted by particle j on particle i is — V<fr(a;j — Xj). 

- When the system is constituted of two elementary particles, in the reference frame attached to the 
center of mass, the dynamics is one-dimensional. The deflection of the particle trajectories from 
straight lines can then be described through explicit formulas (which will be given in Chapter 3). 

- When the system is constituted of three particles or more, the integrability is lost, and in general 
the problem becomes very complicated, as already noted by Poincare [34]. 



1.1. Thermodynamic limit 



In the large N limit, called the thermodynamic limit, individual trajectories become irrelevant, and 
the goal is to describe an average behaviour. 

The Liouville equation relative to the particle system (1.0.1) is 

N N N 

dtfN + ^2 Vi " V **/jv - Va; * ~ x i) ' v «./^ = . 

i=l i=l 3=i 

We use the following notation: for any set of s particles with positions X s := (x\, . . . ,x s ) e R ds and 
velocities V a :— (vi, ...,v s ) £ K ds , we write Zi := (xi,Vi) € R 2d and Z s := (z\, . . . ,z s ) € R 2ds . We 
assume that the probability Jn, referred to as the iV-particle distribution function, satisfies for all 
permutations o of {1, . . . , N}, 

(1-1.1) f N (t,Z„ {N) ) = f N (t,Z N ), 

with Z a ( N ) = (a; CT (i),w CT (i), • ■ • , aV(jv)>?V(iV))- This corresponds to the property that the particles are 
indistinguishable. 
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The average behavior of the particles is then described by the first marginal f N of the distribution 
function /at, defined by 

Z x ) := J f N {t, Z N )dz 2 ...dz N . 

In this framework, in order for the average energy per particle to remain bounded, one has to assume 
that the energy of each pairwise interaction is small. In other words, one has to consider a rescaled 
potential $ e obtained 

- either by scaling the strength of the force, 

- or by scaling the range of potential. 



1.2. Mean field versus collisional dynamics 



According to the scaling chosen, we expect to obtain different asymptotics. 

• In the case of a weak coupling, i.e. when the strength of the individual interaction becomes small 
(of order 1/-/V) but the range remains macroscopic, the convenient scaling in order for the macroscopic 
dynamics to be sensitive to the coupling is: 

N ^ N N 

dtfw + Y Vi ' v *Jn - Y Y ( Xi ~ x j) ■ V Vi /jv = . 

1=1 i=\ 3 = 1 

Then each particle feels the effect of the force field created by all the (other) particles 

1 N f f 

F N (x) = - jj^Vx&ix-Xj) ~ - / / V x $>{x-y)f { x\t,y,v)dydv. 



In particular, the dynamics seems to be stable under small perturbations of the positions or velocities 
of the particles. 

In the thermodynamic limit, we thus get a mean field approximation, that is an equation of the form 

dtf + v- V x f + F ■ V„/ = 
for the first marginal, where the coupling arises only through some average 

F :=-V x $* J fdv. 

An important amount of literature is devoted to such asymptotics, but this is not our purpose here. 
We refer to [10, 38] for pioneering results, to [22] for a recent study and to [19] for a review on that 
topic. 



• The scaling we shall deal with in the present work corresponds to a strong coupling, i.e. to the 
case when the amplitude of the potential remains of size 0(1), but its range becomes small. We shall 
assume throughout this text the following properties for $ (a short-range potential). 



1.2. MEAN FIELD VERSUS COLLISIONAL DYNAMICS 
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Assumption 1.2.1. — The potential $ : R d — > R is a radial, nonnegative, nonincreasing function 
supported in the unit ball ofR d , of class C 2 in {x e R d , < |a:| < 1}. Moreover it is assumed that $ 
is unbounded near zero, goes to zero at \x\ = 1 with bounded derivatives, and that V$ vanishes only 
on \x\ — 1. Finally writing $(x) = <!>( |cc|) we assume that for all p € (0, 1), 

(1.2.1) p$"(p)+2&(p) >0. 

Remark 1.2.2. — We refer to Chapter 3 for a justification of those assumptions, in particular (1.2.1) 
that appears as a sufficient condition to define a scattering cross-section. Condition (1.2.1) can easily 
be checked for a large class of potentials. For instance any potential of the form 3>(p) = p~ k — 1 
for p < 1 is suitable (for k > 1). Potentials smooth at p = 1 can be constructed from that example by 
using a smooth junction ([36]). 



Introduce a small parameter e > corresponding to the typical interaction length of the particles. 
Then in the macroscopic spatial and temporal scales, the Hamiltonian system becomes 

dxi dvi 1 ^ 



and the Liouville equation takes the form 

N N N 



(1.2.3) d t f N + ^v l -V x J N -Y,Y,-V x ^( X ^A -V Vi f N =0. 

i=l i=l 3 = 1 £ V £ / 

With such a scaling, the dynamics is very sensitive to the positions of the particles. 




Figure 1. Instability 

Situations 1 and 2 on Figure 1 are different by a spatial translation of 0(s) only. However in Situation 1, 
particles will interact and be deviated from their free motion, while in Situation 2, they will evolve 
under free flow. 
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1.3. The Boltzmann equation 

Of course, particles move with uniform rectilinear motion as long as they remain at a distance greater 
than e to other particles. In the limit e — » 0, wc thus expect trajectories to be almost polylines. 

Deflections are due to elementary interactions 

- which occur when two particles are at a distance smaller than e, 

- during a time interval of order e (if the relative velocity is not too small), 

- which involve generally only two particles : the probability that a third particle enters a security 
ball of radius e should indeed tend to as e — > in the convenient scaling. We are therefore brought 
back to the case of the two-body system, which is completely integrable (see Chapter 3). 

In order for the interactions to have a macroscopic effect on the dynamics, each particle should undergo 
a finite number of collisions per unit of time. A scaling argument, giving the mean free path in terms 
of N and e, then shows that Ne^ 1 = 0(1). This is the Boltzmann-Grad scaling (see [21]). 

In the limit e ->• with Ne^ 1 = 1, we would like to obtain a kind of homogeneisation result : we 
want to average the motion over the small scales in t and x, and replace the localized interactions by 
pointwise collisions as in the case of hard spheres. 

We shall therefore introduce an artificial boundary (following [27]) so that 
on the exterior domain, the dynamics reduces to free transport, 

- on the interior domain, the dynamics can be integrated in order to compute outwards boundary 
conditions in terms of the incoming flux. Note that such a scattering operator is relevant only if wc 
can guarantee that there is no other particle involved in the interaction. 

The statistical distribution of deflection angles w is then predicted by a function b = b(v — vi,u), the 
collision cross-section, which depends only on the microscopic interaction potential. 

A counting argument leads then to the Boltzmann equation (introduced by Boltzmann in [7]- [8]) : 
' d t f + v-V x f = Q(JJ) 

free transport localized binary collisions 

QifJ) ■■= J J [fiv^fiv'J-f^fiv^biv-vu^dviduj. 

v' , v[ pre-collisional velocities 

The collision term, which acts only on the ^-variable, is constituted of a gain term, corresponding to 
the creation of particles of velocity v by collision between particles of velocities v' and v[ , and of a loss 
term, due to the disappearance of particles of velocity v by collision with particles of velocity v\ . 

Note that the joint probability of having particles of velocity (i/, v[) (respectively of velocities (y, v\)) 
before the collision is assumed to be equal to fit, x, v')f(t, x, v[) (rcsp. to f(t, x, v)f(t, x, Vi)), meaning 
that there is independence. 

1.4. A convergence result 

The goal of this text is to prove the following statement. It is written somewhat loosely, we refer 
to Chapter 7 for a precise statement (see in particular the statement of Theorem 4 page 53). The 
appropriate notion of independence is defined in Section 7.1 and the appropriate notion of convergence 
is defined in Section 7.2. 



1.5. RELATED RESULTS 
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Theorem 1. — Assume that the repulsive potential <& satisfies Assumption 1.2.1. Let fo : R 2d i-> R + 
&e a continuous function, of integral one, and with exponential decay at large energies. Consider the 
system of N particles, initially distributed according to fo and "independent" (in a sense made precise 
in Chapter 7 ) , governed by the system (1.2.2). Then, in the Boltzmann-Grad limit N — > oo, Ne^ 1 <~ 1, 
its distribution function converges to the solution to the Boltzmann equation (1.3.1) with initial data 
fo, in the sense of observables. 



1.5. Related results 

The problem of asking for a rigorous derivation of the Boltzmann equation from the Hamiltonian 
dynamics goes back to Hilbert [24], who suggested to use the Boltzmann equation as an intermediate 
step between the Hamiltonian dynamics and fluid mechanics, and who described this axiomatization 
of physics as a major challenge for mathematicians of the twentieth century. 

We shall not give an exhaustive presentation of the studies that have been carried out on this question 
but indicate some of the fundamental landmarks. First one should mention N. Bogoliubov [5], M. 
Born, and H. S. Green [9], J. G. Kirkwood [28] and J. Yvan [42], who gave their names to the 
BBGKY hierarchy we shall be using extensively in this study. H. Grad was able to obtain in [20] 
a differential equation on the first marginal which after some manipulations converges towards the 
Boltzmann equation. The first mathematical result on this problem goes back to C. Cercignani [11] 
and O. Lanford [31] who proved the propagation of chaos by a careful study of trajectories of a hard 
spheres system, and who exhibited - for the first time - the origin of irreversibility. The proof, even 
though incomplete, is therefore an important breakthrough. The limits of their methods, on which 
we will comment later on - especially regarding the short time of convergence - are still challenging 
questions. 

The argument of O. Lanford was then revisited and completed in several works. Let us mention 
especially the contributions of K. Uchiyama [39], C. Cercignani, R. Illner and M. Pulvirenti [14] and 
H. Spohn [37] who introduced a mathematical formalism, in particular for the existence of solutions 
to the BBGKY hierarchy which turns out to be a theory in the spirit of the Cauchy-Kowalewskaya 
theorem. 

The term-by-term convergence of the hierarchy in the Boltzmann-Grad scaling was studied in more 
details by Cercignani, V. I. Gerasimcnko and D. I. Petrina [13] : it indeed requires refined estimates 
on the set of "pathological trajectories" , i.e. trajectories for which the Boltzmann equation does not 
provide a good approximation of the dynamics. 

The method of proof was extended 

- to the case when the initial distribution is close to vacuum, in which case global in time results may 
be proved [14, 25, 26]; 

- to the case when interactions are localized but not pointwise [27]. Because multiple collisions are 
no longer negligible, this requires a careful study of clusters of particles. 

Many review papers deal with those different results, see [17, 35, 41] for instance. Our goal here is 
to provide an elementary and self-contained presentation, which includes all the details of the proofs, 
especially concerning convergence which to our knowledge is not completely written anywhere, and 
concerning reasonable assumptions that can be made on the potential. 
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CHAPTER 2 



FORMAL DERIVATION OF THE BOLTZMANN EQUATION FOR 

HARD-SPHERES 



This chapter is intended to familiarize the reader with the methods and notation related to BBGKY 
and Boltzmann hierarchies. We present formally, in the situation of hard-spheres, the passage from the 
Liouville equation associated with the TV-particle flow to the BBGKY hierarchy, and then the limit to 
the Boltzmann equation. All the results stated here will be proved - in the more complicated case of 
nonlocal interactions via a potential - in the next chapters. Readers already familiar with the subject 
may skip this chapter altogether. 



2.1. The TV-particle flow 

We consider TV particles, the motion of which is described by N positions (xi, . . . , Xn) and N ve- 
locities («!,..., vn), each in R d . Denoting by Z N := {z Xl . . . , z N ) the set of particles, each parti- 
cle Zi := (xi, Vi) € R 2d is submitted to free flow on the domain 



velocities are precollisional). 

Contrary to the potential case studied in the next chapters, it is not obvious to check that (2.1.1) defines 
a global dynamics, at least for almost all initial data. Note that this is not a simple consequence of 
the Cauchy-Lipschitz theorem since the boundary condition is not smooth, and even not defined for 
all configurations. In the presence of a potential, we shall prove in this text that the set of trajectories 
involving multiple collisions has zero measure. Let us prove this result for the hard sphere dynamics: 
we call pathological a trajectory such that 




and bounces off the boundary &D N according to the laws of elastic reflection: 




v ln ) < (meaning that the ingoing 



- either there exists a collision involving more than two particles, hence the boundary condition is not 
well defined; 
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- or there are an infinite number of collisions in finite time so the dynamics cannot be globally defined. 

In [2], Proposition 4.3 it is stated that outside a negligible set of initial data there are no pathological 
trajectories; the complete proof is provided in [1]. Actually the setting of [1] is more complicated than 
ours since the case of an infinite number of particles is considered. The arguments of [1] can however 
be easily adapted to our case to yield the following result, whose proof we detail for the convenience 
of the reader. 

Proposition 2.1.1. — Let N,e be fixed. The set of initial configurations leading to a pathological 
trajectory is of measure zero in R 2dN . 

We first prove the following elementary lemma. 

Lemma 2.1.1. — Let p,R> be given, and 8 < e/2. Define 

I := {Zn £ x Bft I one particle will collide with two others on the time interval [0, <5] j . 
Then \I\ < C(N,e,R)p d( ^ N ^5 2 . 

Proof. — Just notice that L is embeded in 

^Z N e Bp x Br /3{i,j,k} distinct, \xi - Xj\ e [s, s + 2RS] and \xi - Xk\ € [e, e + 2R5] j , 
and the lemma follows directly. □ 

Proof of Proposition 2.1.1. — Now let R > be given and fix some time t > 0. Let S < e/2 be a 
parameter such that t/S is an integer. 

Lemma 2.1.1 implies that there is a subset I (6, R) oiB^xB^ of measure at most C(N,e, R)R d( ~ N ~ 2) 5 2 
such that any initial configuration belonging to (B R x B R ) \ I (S, R) generates a solution on [0, S] such 
that each particle encounters at most one other particle on [0,(5]. 

Now let us start again at time 5. We recall that in the velocity variables, the ball of radius R in R dN 

is stable by the flow, whereas the positions at time S lie in the ball B^ +RS . Let us apply Lemma 2.1.1 

again to that new initial configuration space. Since the measure is invariant by the flow, we can 

construct a subset 1^5, R) of the initial positions B% x^, of size C(N, e, R)R d ( N ~V (1 + 8) d ( N -^5 2 

such that outside Iq U 1\ {8, R) , the flow starting from any initial point in B R x B R is such that each 

particle encounters at most one other particle on [0, 5], and then at most one other particle on [6, 25]. 

t/5-i 

We repeat the procedure t/S times: we construct a subset Ig(t,R) := Ij(S,R) of B R x B R , of 

3=0 

measure 

t/S-i 

\I s (t,R)\<C(N,s,R)R d ^ N - 2) 5 2 {l+jS) d(N - 2) <C(N,R,t,e)S 

3=0 

such that for any initial configuration in B R x B R outside that set, the flow is well-defined up to time t. 
The intersection I(t,R) := O Is(t,R) is of measure zero, and any initial configuration in B R x B R 

8>0 

outside I(t, R) generates a well-defined flow until time t. Finally we consider the countable union of 
those zero measure sets / := [J I(t n , R n ) where t n and R n go to infinity, and any initial configuration 

n 

in H 2dN outside / generates a globally defined flow. The proposition is proved. □ 



2.2. THE LIOUVILLE EQUATION AND THE BBGKY HIERARCHY 
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2.2. The Liouville equation and the BBGKY hierarchy 

The Liouville equation relative to the particle system (2.1.1) is 

N 

(2.2.1) d t f N + ^Vi-V x J N = Q in V t 



N 
»=1 

with the boundary condition /jv(t, ZJJ 1 ) = /jv(t, ■Zjv**)- We assume from now on that /jv is invariant 
by permutation in the sense of (1.1.1), meaning that the particles are indistinguishable. 



One can associate with this Liouville equation a hierarchy of equations, satisfied by the marginals 
/«'(^s) / /iv(i, Z s , z s+ i, . . . , z n )1 Zn ev n dz s+ i . . . dz N . 

jR!«i(»->) 

Let us derive this hierarchy formally. We integrate (2.2.1) over the (N — s) last variables, and we first 
notice that 

/ 9 t fN{t, Z s ,z s+ i, . . . , z N )lz N ev N dz s+ i . . . dz N = d t f^\t, Z s ) . 

J R 2d(JV- s ) 

Next we compute 

N f 

V" / Vi ■ V x J N {t, Z N )lz N ev N dz s+ i . . .dz N 

using Green's formula. The boundary terms involve configurations with at least one pair satisfy- 
ing 1 < i < N and s + 1 < j < N, with \xi — Xj \ = e. According to the previous section we may neglect 
configurations where more than two particles collide at the same time, so the boundary condition is 
well defined. For any i e {1, . . . , N} and any j e {s + 1, . . . , N}, we denote by the outwards 

normal at any point of the boundary and we define 

Z N (i,j) := [Z N € K 2dN , \xi - Xj \ = e} . 

X ' — X ' 

Recalling that v' l 'i := — we then obtain, using the invariance of fn by permutation, 

| X i X j | 

N 



V" / Vi - V x J N (t, Z N )lz N ev N dz s+1 . . . dz N 

~[ J-R2d{N-s) 



i=l j=li=s+1 ^(i,j) 
with da 1 ^ the surface measure on Ejv(«,j), induced by the Lebesgue measure. 
By symmetry this gives 

AT 

Vi ■ V Xi /jv(t, Z N )tz N ev N dz s+ i . . . dzn 



' N (t,Z N ) da% 3 , 



2d(N-s) 



= -E« i -v X4 /W(t,z.) + ^£/ 

~[ V2 i=1 Jx N (i,s+i) 

It remains to define the collision operator 



v l > s+1 -(v^-Vi) f N {t,Z N )d<j^ s - 



(2.2.2) C s . s+1 f^ +1 \t,Z s+1 ) := (N - S )f" [ v^ 1 ■ (v s+1 - v { ) f% +1 \t, Z s , z s+1 ) dadv s+1 

l=1 JS e (x t )xTL<i 
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where S e (xi) is the sphere of radius e centered at Xi and da is the surface measure on that sphere and 
in the end we obtain the BBGKY hierarchy 

(2.2.3) + v i ■ V -J ( n = C s,s+if { N +1) in V s , 

l<i<s 

with the boundary conditions (2.1.1). 



2.3. The Boltzmann hierarchy and the Boltzmann equation 

Starting from (2.2.3) we now consider the limit N — > oo under the Boltzmann-Grad scaling Ne^ 1 = 1. 
The Duhamel formulation for (2.2.3) writes 

/&>(*) = D s (t).4 s) + fT> s {t-T)C s , s+1 f^ +1 \r)dT, 
Jo 

where D s (i) denotes the s-particle flow on V s with the boundary conditions (2.1.1). 

Because of the scaling assumption, the collision term C s ^ + if^ s+1 \Z s ) is approximately equal to 

s . 

-{N - s)e d ~ 1 V] / v- (v s+1 -Vi)f^ +1) (Z s ,Xi+ su,v s+ i )dvdv s+ i 

which we may split into two terms, depending on the sign of v ■ (v s+ i — Vi): 

s „ 

y2 / (v-(Vs+i-Vi)) fx +1 \z s ,Xi +ev,v s+1 )dvdv s+1 

s „ 

(v-(v a +i-Vi)) f^ +1 \z s ,Xi + eu,v s+ i)dudv s+ i. 

~^Jsf- 1 xR dV 

Recall that pre-collisional particles are particles (xi,v£) and (x s +i,v s +i) for which 

(x s+1 - x^ ■ (v s+ i - < . 

The case when (ar s +i — Xi) ■ (u s +i — u») > is called the post-collisional case. Consider a set of 
particles Z s+ i such that (xi,Vi) and (a; s +i,w s +i) are post-collisional. Provided that there is no other 
particle among these s + 1 which has undergone a collision on a short time interval (which is almost 
sure in the limit e — > 0), we have 

f { N +1 \t,Z s ,x s+1 ,v s+1 ) = f^ +1) (t,Z*,x* +1 ,v* +1 ) 

where (z* , z* s+l ) is the pre-image of (zi, z s +i) by (2.1.1). Then neglecting the small spatial translations 
in the arguments of f^ +1 \ we obtain the following asymptotic expression for the collision operator at 
the limit: 

Cs,s+J {s+1 \t, Z s ) :=J2[ l Hvs+1 _ Vt)>0 v ■ (v s+1 - Vi ) 
»=i J 

x (/ (s+1) (Mi,i>i,-- -,Xi,v*,.. .,x s ,v s ,Xi,v* +1 ) - / (s+1) (i, Z s ,Xi,v a+1 fjdudv a+1 . 
The asymptotic dynamics are therefore governed by the following Boltzmann hierarchy: 

(2.3.1) /W(t) = S s (t)ti s) + rS s (i-r)CL +1 / (s+1) (T)dr. 

Jo 

where S s (t) denotes free-flow. 



2.3. THE BOLTZMANN HIERARCHY AND THE BOLTZMANN EQUATION 
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s 

Note that if f( s \t,Z s ) = J^[/(A z i) (meaning f^ s \i) is tensorized) then / satisfies the Boltzmann 

»=i 

equation (1.3.1), where the cross-section b is simply b(v s+ i — v i7 uj) :— l u .(„ s+1 _„ i )>o w • (w s +i — ^j). 



CHAPTER 3 



TWO-PARTICLE INTERACTIONS 



In the case when the microscopic interaction between particles is governed by a short-range repulsive 
potential, collisions are no more instantaneous and pointwise, and they possibly involve more than two 
particles. Our analysis in Chapters 8 to 10 shows however that the low density limit Ne^ 1 — > requires 
only a description of two-particle interactions, at the exclusion of more complicated interactions. 

In this chapter we therefore study precisely, following the lines of [12], the Hamiltonian system (1.2.2) 
for N — 2. The study of the reduced motion is carried out in Section 3.1, while the scattering map is 
introduced in Section 3.2, and the cross-section, which will play in important role in the Boltzmann 
hiearchy, is described in Section 3.3. 



3.1. Reduced motion 

Wc first define a notion of pre- and post-collisional particles, by analogy with the dynamics of hard 
spheres: 

Definition 3.1.1. — Two particles zi 7 z 2 are said to be pre-collisional if they belong to the artificial 
boundary and their distance is decreasing: 

\xi - x 2 \ =£, (vi - v 2 ) ■ (xi - x 2 ) < 0. 

Two particles z\,z 2 are said to be post-collisional if they belong to the artificial boundary and their 
distance is increasing: 

\xi -x 2 \=s, Oi - v 2 ) ■ (x 1 - x 2 ) > 0. 

We consider here only two-particle systems, and show in Lemma 3.1.2 that, if z\ and z 2 are pre- 
collisional at time i_ , then there exists a post-collisional configuration z[ , z 2 , attained at t + > t-. 
Since V<f>(x) vanishes on {|x| > e}, the particles z\ and z 2 travel at constant velocities v[ and v 2 for 
ulterior (t > t + ) times. 

Momentarily changing back the macroscopic scales of (1.2.2) to the microscopic scales of (1.0.1) by 
defining r := (t — t-)/e and y(r) := x/s(t), w(t) — v(t), we find that the two-particle dynamics is 
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governed by the equations 
(3.1.1) 




(3.1.2) ^-(w 1 +w 2 )=0, 4~ {-Awi +w 2 ) 2 + hw! -w 2 ) 2 + <P(yi - y 2 . 



whence the conservations 

Tr {Wl+ '-' ~> dr\^ 1 4 

From (3.1.2) we also deduce that the center of mass has a uniform, rectilinear motion: 

(3.1.3) (yi + ife)(r) = (?a + y 2 ){0) + t( Wi + w 2 ) , 

and that pre- and post-collisional velocities are linked by the classical relations 

(3.1.4) w[ + w' 2 = wi + w 2 , Kl 2 + KI 2 = HI 2 + M 2 - 

A consequence of (3.1.1) is that (5y,Sw) := (t/i — y 2 ,wi — w 2 ) solves 

(3.1.5) ^5y = Sw, ^-5w = -2V$(5y) . 

dr dr 

We notice that, $ being radial, there holds 

^-(5y A Sw) — Sw A Sw - 28y A V$(<5?y) = , 
dr 

implying that, if the initial angular momentum Syo A Swq is non-zero, then 5y remains for all times in 
the hyperplane orthogonal to Syo A Swo- In this hyperplane, introducing spherical coordinates (p,f) 
in R+ x S^" 2 , such that 

Sy — pe p and Sw = pe p + pipe v 
the conservations of energy and angular momentum take the form 

1 -(pi + (p<pf) + 2<l>(p)= 1 -\5w \\ 

p 2 M = \$vo a Sw \ , 

implying p > for all times, and 

(3.1.6) p 2 + V(p,£ ,J )=£o, * : =^<L + 4<I>(p), 
where we have defined 

(3.1.7) £q := \5w n \ 2 and Jo := \Syo A <5wo|/|^wo| = sina , 

which are respectively (twice) the energy and the impact parameter, ir — a being the angle between Swq 
and Syo (notice that a > tt/2 for pre-collisional situations). In the limit case when a = 0, the movement 
is confined to a line since ip = 0. 

We consider the sets corresponding to pre- and post-collisional configurations: 

(3.1.8) S ± := {(Sy,6w) E S^ 1 x R d , ±Sy ■ Sw > 0} , 

where S^ 1 is the unit sphere centered at the origin in R d ; in spherical coordinates pre-collisional 
configurations correspond to p = 1 and p < while post-collisional configurations are correspond to 
p = 1 and p > 0. 



3.1. REDUCED MOTION 



15 



Lemma 3.1.2 (Description of the reduced motion). — For the differential equation (3.1.5) 
with pre-collisional datum (Syo,dwo) € S~ , there holds \6y(T)\ > p* for all r > 0, with the notation 

(3.1.9) p* =p*(£ Q ,Jo) := max {p € (0,1), *(p,f 0) >7o) = £o} , 
and /or defined by 

(3.1.10) n:=2/ (£ -*(p,£ ,Jo)r 1/2 ^, 
£/ie configuration is post-collisional (p = 1, p > Oj at t = t*. 

Proof. — Solutions to (3.1.6) satisfy p = t(p)(£ — ^(p)) 1 ^ with i(p) = ±1, possibly changing 
values only on {vp = £o}, by Darboux's theorem (a derivative function satisfies the intermediate value 
theorem). The initial configuration being pre-collisional, there holds initially i = — 1, corresponding 
to a decreasing radius. The existence of p* satisfying (3.1.9) is then easily checked: we have |o"yo| = 1 
and 8yo ■ Swq ^ 0, so there holds \I/(1, £q, J$) < £q, and 'F is increasing as p is decreasing. The set 
{r > 0,p(r) > p*} is closed by continuity. It is also open: since <F is nonincreasing, then d p ^ ^ 
everywhere and in particular at (p*, £g, Jo). So £q — 'F changes sign at p*, which forces, by (3.1.6), 
the sign function t to jump from — to + as p reaches the value p* from above. This proves p > p* by 
connexity. The minimal radius p = p* is attained at t*/2, where is defined by (3.1.10), the integral 
being finite since d p ^ does not vanish. Assume finally that for all r > 0, there holds p(r) < 1. Then 
on [r + /2, +oo), p is increasing and bounded, hence converges to a limit radius, which contradicts the 
definition of p*. This proves p = 1 at t = r*, a time at which p > 0, since i has jumped exactly once, 
by definition of p* . □ 

Remark 3.1.3. — Denoting A : (y,w) — > (y,—w), and <f> t : R 2d — > R 2d the flow of (3.1.5), we find 
that <j)-t = A o cj) t o A, implying A o <j) t o A o (f> t = Id, and time-reversibility of the two-particle dynamics. 




Figure 2. Reduced dynamics 

The reduced dynamics is pictured on Figure 2, where the half-deflection angle 9 is the integral of the 
angle <p as a function of p over [p* , 1] : 

(3.1.11) 6= f ^^(So-^(p,£ ,l )r 1/2 dp, 

J P , P 
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With the initialization choice ipo = 0, the post-collisional configuration is (p, <p)(t*) — (1,29); it can 
be deduced from the pre-collisional configuration by symmetry with respect to the apse line, which by 
definition is the line through the origin and the point of closest approach (5y(r*/2), 5w(t*/2)). The 
direction of this line is denoted to G S 



3.2. Scattering map 

We shall now define a microscopic scattering map <t that sends pre- to post-collisional configurations: 

<7 : (6yo,8w ) G S~ — > (8y(T*),6w(T*)) = 4> Tt {5y ,8w ) G S + . 

By uniqueness of the trajectory of (3.1.5) issued from (8yo,5wo) (a consequence of the regularity 
assumption on the potential, via the Cauchy-Lipschitz theorem), the scattering is one-to-one. It is also 
onto, by Remark 3.1.3: the pre-image of (Sy, Sw) G S + by the scattering is I o cf> Tt (5y, —Sw) G S~ . 

Back in the macroscopic variables, we now define a corresponding scattering operator for the two- 
particle dynamics. In this view, we introduce the sets 

St := {(z u z 2 ) G R 4d , \x!-x 2 \ =e, ±(an - x 2 ) ■ («i - «a) > o} . 

We define, as in (3.1.7), 

/noil CI 1 2 AT |(X1 — X 2 ) A (Vl — ^ 2) I 

(3.2.1) £ = \vi — v 2 \ and J := ; ; =: sin a. 

e\vi - v 2 \ 

Definition 3.2.1 (Scattering operator). — The scattering operator is defined as 

a e : {x 1 ,v 1 ,x 2 ,v 2 ) G S~ — > (x'^v'^x'^v^) G «S+ , 

where 



(3.2.2) 



where is the microscopic interaction time, as defined in Lemma 3.1.2, (5y(r*), 6w(t*)) is the micro- 
scopic post-collisional configuration: ((h/(r*), 5w(t*)) = &o(( x i — x 2 )/s, v\ — v 2 ), and cu is the direction 
of the apse line. Denoting by v := (x\ — x 2 )/\xi — x 2 \ we also define 

o-o(v,v\,v 2 ) := (v',v[,v' 2 ) . 

The above description of {x\,v' x ) and (x' 2 ,v 2 ) in terms of to is deduced from the identities 

Sv(t*) = 5v - 2ui ■ Sv u> and Sy(n) = -Sy + 2ui ■ Sy co 
in the reduced microscopic coordinates. 

By d p ^ 7^ in (0, 1) and the implicit function theorem, the map [£,J) — > p*(£,J) is C 2 just like \E". 
Similarly, G C 2 . By Definition 3.2.1 and C 1 regularity of V<I> (Assumption 1.2.1), this implies that 
the scattering operator a e is C 1 , just like the flow map of the two-particle scattering (denoted </> m 



x\ 


■ = ^{xi + x 2 ) + + v 2 ) + E -8y{T,\ 


) = -xi + LO ■ 


(Xl 


- x 2 )oj + ^-{vi + v 2 ) , 


x 2 


■= ^{xi + x 2 ) + ^j-{vi + v 2 ) - |(5y(r*; 


) = -x 2 - uj ■ 


(xi 


- x 2 )oj + ^-{vi + v 2 ) , 


v[ 


:= i(ui + v 2 ) + ^5w(t„) =vi-w-(vi 


-v 2 )u>, 






v' 2 


■ = ^{vi +v 2 ) - ^5w(t*) = v 2 +lu ■ (vx 


-v 2 )uj, 







3.2. SCATTERING MAP 
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Remark 3.1.3). The scattering a £ is also bijective, for the same reason that the microsopic scattering 
<T is bijective; the inverse map is a^ 1 := / o <j £ o I , with notation introduced in Remark 3.1.3. 

For any s £ N*, R > 0, we denote B S R := {V s £ K ds , \V S \ < R} where | • | is the euclidean norm; we 
often write Br := B R . 

Proposition 3.2.1. — Let R > be given and consider 

S t,R : = {( z i' z 2) G (R d x B R f, \x!-x 2 \ = e,± (v 1 -v 2 ) ■ (x 1 -x 2 ) > o} . 
The scattering operator a e is a bijection from S~ R to S^~ R . 

The macroscopic time of interaction % = %(£q, Jo) '■= £T*, where is defined in (3.1.10), is uniformly 
bounded on compact sets of R + \ {0} x [0, 1]. 

Proof. — We already know that a £ is a bijection from S~ to <S+. By (3.1.4), it also preserves the 
velocity bound. Hence a E is bijective S~ R — > S^~ R . 

Now given £ > and Jo € [0, 1], we shall show that r* can be bounded by a constant depending only 
on £ . Since $(p*) < £o/4, then p* > <J> _1 (f /4). Let us then define i £ (0, 1) by 

so that p 2 > 8if t . 

On the one hand it is easy to see, after a change of variable in the integral, using 

Us, J ,p)) = 2£ °f- W(p) > 2£ ^l > 2£ J 2 , 

dp p A p 6 

that there holds the bound 



1 f 

£oJq Jo 



dy < 2^i~j2 



Vv ~ Jo 2 V£o~ 

So if Jo > *0) we find that 

2 16 

T* < 



On the other hand for Jo < io we define 7 := $~ 1 (£o/8) and we cut the integral defining r* into two 
parts: 

r. = ri 1} + ri 2) with t* (1) = 2 T (£„ - Jo, p))~ 1/2 ^ • 

Notice that since p 2 > 8ig and J < i , then £ /4 - EoJol^pl > 7f /32 > £ /8 so 

».=*-( T -Tf)^-'(!H' 

The first integral t* is estimated using the fact that $' does not vanish outside 1 as stated in 
Assumption 1.2.1: defining 

we find that on [10,7], 



M(*):= inf |$'(p)|>0, 

*o<P<7 



d IF 7 2 

-f (£„ - *(£o, Jo,p)) = -^p- - 4$'(p) > 4M($) 
dp p 01 
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SO 

Tl U,2\ 2 



(1) (fo/2 - goJ-pW) 2 < V^O 
T * - M($) " 72M($) 

For the second integral we estimate simply 

(2) 2 2 4^2 

' * _^ 1 _^ 



(£ /2-foJo 2 /7 2 ) 5 (£ /2-fo/8) 5 
The result follows. □ 

Remark 3.2.2. — //$ is convex then M(<1>) = |$'(7)|. Moreover if $ is of the type exp(— - - 2 ) 
then the proof of Proposition 3.2.1 shows that r* may be bounded from above by a constant of the order 
of Cj ^(1 + log e ) i/ £ o > e . 



3.3. Scattering cross-section and the Boltzmann collision operator 

The scattering operator in Definition 3.2. f is parametrized by the impact parameter and the two 
ingoing (or outgoing) velocities. However in the Boltzmann limit the impact parameter makes no 
longer sense: the observed quantity is the deflection angle or scattering angle, defined as the angle 
between ingoing and outgoing relative velocities. The next paragraph defines that angle and as well 
as the scattering cross-section, and the following paragraph defines the Boltzmann collision operators 
using that formulation. 



3.3.1. Scattering cross-section. — With notation from the previous paragraphs, the deflection 
angle is equal to n — 29 where 6 := a + 9, the angle a being defined in (3.2.1) and 9 being defined in 
cqrefthcta, so that 



6 = Q(£ ,Jo) := arcsinJo + Jo 




The following result, and its proof, are due to [36]: 

Lemma 3.3.1. — Under Assumption 1.2.1 and for all £ > 0, the function Jo i->- 0(£ o , Jo) € [0, 7r/2] 
satisfies Q(£o,0) = and is strictly monotonic: dj @ > for all J e (0,1). Moreover, it satis- 
fies lim dj & G (0, oo ] and lim dj Q = 0. 

Jo— >0 " Ja— >1 



Proof. — An energy £ > being fixed, the limiting values @(£o, 0) = and Q(£o, 1) = tt/2 arc found 
by direct computation. To prove monotonicity, the main idea of Saffirio and Simonclla is to use the 
change of variable 

sm 2 if := — ^ + -f 
£q P 2 

which yields 

/2 sin ip 

resin Jo — ~ £oJq 



3.3. SCATTERING CROSS-SECTION AND THE BOLTZMANN COLLISION OPERATOR 
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Computing the derivative of this expression with respect to j7o gives 



dJo yw ""' TW^V 

J arcsin J a 

where ip is defined by 



p<S>"(p) + 2<£'(p) + -^(<S>'(p)) 2 )<bp 



P 2 to 

In view of the formula giving dj Q, it turns out assumption (1.2.1) implies dj Q > for all Jq £ (0, 1), 
and also the limits 

lim dj Q £ (0,oo] and lim 9^0 = 

Jo— >0 J'o— >1 

as soon as $'(1) = (if not then lim SjqO = oo). The result follows. □ 

c7b — >•! 

Remark 3.3.2. — Note that one can construct examples that violate assumption (1.2.1) and for which 
monotonicity fails, regardless of convexity properties of the potential <I> ([36]). 

By Lemma 3.3.1, for each £ we can locally invert the map 0(£oj •)> and thus define J7b as a smooth 
function of £q & n d 0- This enables us to define a scattering cross-section (or collision kernel), as follows. 




Figure 3. Spherical coordinates 

For fixed x\, we denote do\ the surface measure on the sphere {y £ R d , \y — xi\ = e}, to which x 2 
belongs. We can parametrize the sphere by (a,ip), with ip £ S d ~ 2 , where a is the angle defined in 
(3.2.1). There holds 

da i = £ d ~ x (sin a) d ~ 2 dudi\) . 
The direction of the apse line is to — (Q,ip), so that, denoting duj the surface measure on the unit 
sphere, there holds 

(3.3.1) doj = (sin e) d ~ 2 d@dip. 

By definition of a in (3.2.1), there holds 

(xi - x 2 ) ■ (vi - v 2 ) = e\vi - v 2 \ cos a, 
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SO that 

- (xi — x 2 ) • (vi — v 2 ) dcr\ = e d ~ 1 \vi — v 2 \ cos a (sma) d ~ 2 dadip 

= s d - 1 \v 1 -v 2 \jt 2 dJodip, 
where in the second equality we used the definition of J in (3.2.1). This gives 

(3.3.2) l - ( Xl - x 2 ) ■ { Vl - v 2 ) da x = e d ->! - v 2 \J d - 2 d@J dOdiP , 
wherever d® Jo is defined, that is, according to Lemma 3.3.1, for J e [0, 1). 

Definition 3.3.3. — The scattering cross-section is defined for \v\ — v 2 \ > and O G (0, 7r/2] by 

(3.3.3) b(\vi-V2\,Q) := \v 1 -v 2 \J d - 2 deJo(^me) 2 - d . 
Abusing notation we shall write b(\vi — v 2 \,Q) — b(\vi — v 2 \,uj). 

By Lemma 3.3.1, the cross-section b is a locally bounded function of the relative velocities and scattering 
angle. 

3.3.2. Scattering cross-section. — The relevance of b is made clear in the derivation of the Boltz- 
mann hierarchy, where we shall use the identity 

(3.3.4) * (x 1 - x 2 ) ■ (vi - v 2 ) da 1 = e d ~ 1 b{\v 1 - v 2 \,w)dco , 

derived from (3.3.1), (3.3.2) and Definition 3.3.3. As in Chapter 2 (see in particular Paragraph 2.3), 
we can formally derive the Boltzmann collision operators using this formulation: we thus define 

i=i J 

(3.3.5) x ^/ (s+1) (t,a;i,wi, ...,Xi,v*,.. . , x s , v s , x h v* +1 ) - / (s+1) (i, Z s , x u v s+ i)^dvdv s+1 , 
where {v*,v* +1 ) is obtained from (vi,v s+ i) by applying the inverse scattering operator a^ 1 , using 



_ w x i x s+l \ f %i X s +1 * * \ 



z s +i| 

This can also be written using the cross-section: 

C 5 ° s+ i/ (s+1) (*,^):=E / Kh-«2|.w) 

x (f (s+1) (t,x 1 ,v 1 , ...,Xi,v*,.. .,x s ,v s ,Xi,v* +1 ) - f {s+1) (t,Z s ,Xi,v s+1 )jdujdv s+1 . 

Remark 3.3.4- — It * s n °t possible to define an integrable cross-section if the potential is not com- 
pactly supported, no matter how fast it might be decaying. This issue is related to the occurrence of 
grazing collisions and discussed in particular in [41], Chapter 1, Section 1.4. However it is still possible 
to study the limit towards the Boltzmann equation, if one is ready to change the formulation of the 
Boltzmann equation by renouncing to the cross-section formulation ([36]). 

The question of the convergence to Boltzmann in the case of long-range potentials is a challenging open 
problem; it was considered by Desvillettes and Pulvirenti in [15] in the linear case, while Desvillettes 
and Ricci studied grazing collisions in [16]. 



CHAPTER 4 



THE BBGKY HIERARCHY 



The main goal of this text is to extend the formal strategy described in Chapter 2 for hard spheres to 
general short-range potentials, then to rigorously justify all the steps of the convergence proof. This 
necessitates the derivation of bounds for Jn that do not depend on N. 

Our starting point is the Liouville equation (1.2.3) satisfied by the A^-particle distribution function /jy. 
We reproduce here equation (1.2.3): 

(4.0.1) d t f N + £ Vi-V Xi f N - ]T ^(°^A-V v J N = 0. 

l<i<N l<i^j<N ^ ' 

The arguments of /jv in (4.0.1) are (t, Zn) € R + x fijy, where 

£l N := [Z N e H 2dN ,Vi^j,x t ^x 3 }. 

The classical strategy to obtain a kinetic equation such as (1.3.1) is to write the evolution equation for 
the first marginal of the distribution function f N , namely 

f$(t,zi):= / f N (t,z 1 ,z 2 ,...,z N )dz 2 ...dz N . 

J R 2d(JV-l) 

The point to be noted is that the evolution of fff depends actually on fffl because of the quadratic 

(2) 

interaction imposed by the force F = — V<i>. And in the same way, the equation on f N ' depends 
on fffi . Instead of a kinetic equation, we therefore obtain a hierarchy of equations involving all the 
marginals of /at 

(4.0.2) f$(t,Z a ):= f N (t,Z s ,z s+1 ,...,z N )dz s+1 ---dz N . 

J R 2d(N-s) 

In Section 4.1 it is shown that due to the presence of the potential, and contrary to the hard spheres 
case described in Chapter 2, it is necessary to truncate those marginals away from the set Hn- An 
equation for the truncated marginals is derived in weak form in Section 4.2. In order to introduce 
adequate collision operators, the notion of cluster is introduced and described in Section 4.3. Then 
collision operators are introduced in Section 4.4, and finally the integral formulation of the equation 
is written in Section 4.5. 
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4.1. Truncated marginals 

From (4.0.1), we deduce by integration that the untruncated marginals defined in (4.0.2) solve 
d t .f ( °\t, Z„) + • V x J%\t, Z s )--j^F f x i^ x l) . V vt f%\t, Z s ) 



e 

1=1 i,j = l 

(4.1.1) '** 



N - 



E/ F ( - £ XS+1 ) -V Vi ti +1 \t,Z s ,z s+1 )dz s+1 . 



There are several differences between (4.1.1) and the BBGKY hierarchy for hard spheres (2.2.2)-(2.2.3). 
One is that the transport operator in the left-hand side of (4.1.1) involves a force term. Another is that 
the integral term in the right-hand side of (4.1.1) involves velocity derivatives. Also, that integral term 
is a linear integral operator acting on higher-order marginals, just like (2.2.2), but, contrary to (2.2.2), 
is not spatially localized, in the sense that the integral in x s+ i is over the whole ball B(xi,e), as 
opposed to an integral over a sphere in (2.2.2). 

This leads us to distinguish spatial configurations in which interactions do take place from spatial 
configurations in which particles are pairwise at a distance greater than e, by truncating off the 
interaction domain {Z N , \ Xi — Xj\ < e for some i ^ j} in the integrals defining the marginals. For the 
resulting truncated marginals, collision operators will appear as integrals over a piece of the boundary 
of the interaction domain, just like in the case of hard spheres. The scattering operator of Chapter 3 
(Section 3.2) will then play the role that the boundary condition plays in the case of hard spheres, as 
sketched in Chapter 2. 

Suitable quantities to be studied are therefore not the marginals defined in (4.0.2) but rather the 
truncated marginals 

(4-1.2) f^\t,Z s ):= f N (t,Z s ,z s+1 ,...,z N ) 1] a, \ >e dz s+1 ■ ■ ■ dz N , 

ie{s+i,...,w} 

where | • | denotes the euclidean norm. Notice that 

(7n ) ~ /^XMi) = f N (t, Zl ,z 2 ,...,z N ) | I (1 - TL\ Xl _ x \ >e ) dz 2 ■ ■ ■ dz N 

so that 

(4-1.3) ||(7£> - f^)(t)\\ LX{R ^ < C(N - l)s d \\f^(t)\\ L ^ {n2) . 

We therefore expect both functions to have the same asymptotic behaviour in the Boltzmann-Grad 
limit Ne^ 1 = 1. This is indeed proved in Lemma 7.1.3 and Corollary 7.2.3 in Chapter 7. 

Given 1 < i < j < N, we denote dZu^, the 2d(j — « + l)-dimensional Lebesgue measure dzidzi+i . . . dzj, 
and dX^j) the d(j — i + l)-dimensional Lebesgue measure dXidXi + \ . . . dxj. We also define 

(4.1.4) V% := [X N G R dN , € [l,s] x [s + 1,7V], \x t - Xj \ > e} , 

where [1, s] is short for [1, s] fl N = {k £ N, 1 < k < s}. Then the truncated marginals (4.1.2) may be 
formulated as follows: 

(4-1.5) f$(t,Z s )= / f N (t,Z s ,z s+ i,...,z N ) j| 1 1 a, x \ >e 1x N eD% dZ s+ i N . 

je{s+i,...,N} 



4.2. WEAK FORMULATION OF LIOUVILLE'S EQUATION 
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The key in introducing the truncated marginals (4.1.5), following King [27], is that it allows for a 
derivation of a hierarchy that is similar to the case of hard spheres. The main drawback is that 
truncated marginals arc not actual marginals, in the sense that 

(4.1.6) 7 ( n\z s )^ f l B (X s+1 )f^ +1 \z s ,z s+1 )dz s+1 , 

for any B C R d( - S+1 \ in particular if B = R d ( s+1) , simply because V S N is not included in T)^ 1 . 
Indeed, conditions \xj — x s+ i\ > £, for j < s, hold for X N € V S N , but not necessarily for X N € . 
Furthermore, V S N intersects all the T> s ^~ m , for m e [1, N — s]. A consequence is the existence of higher- 
order interactions between truncated marginals, as seen below in (4.4.8). Proposition 5.3.1 in Chapter 5 
states however that these higher-order interactions are negligible in the Boltzmann-Grad limit. 



4.2. Weak formulation of Liouville's equation 

Our goal in this section is to find the weak formulation of the system of equations satisfied by the 
family of truncated marginals (//v^) se [i N ^ defined above in (4.1.5). From now on we assume that /jv 
decays at infinity in the velocity variable. 

Given a smooth, compactly supported function <f> defined on R + x R 2ds and satisfying the symmetry 
assumption (1.1.1), we have 

, , / (d t fN + f^v i -V x J N --f^Y j F( X ±^\.V v J N ){t,Z N ) 

(4.2.1) 7r + xr— V ^ \ e ) J 

x <f>(t, Z s )lx N ev° N dZ N dt = . 

Note that in the above double sum in % and j, all the terms vanish except when € [l,s] 2 and 
when G [s + l,iV] 2 , by assumption on the support of F. 

We now use integrations by parts to derive from (4.2.1) the weak form of the equation in the 
marginals /j^ . On the one hand an integration by parts in the time variable gives 

/ d t f N (t, Z N )<f>(t, Z s )l XNeV s dZ N dt = - [ f N (0, Z N )<t>(0, Z s )\ XNeV s dZ N 

/ f N (t,Z N )d t <j)(t,Z s )lx N ev%dZ N dt, 

hence, by definition of 

[ d t f N (t, Z N )<f>(t, Z s )l XN ev° dZ N dt = - [ /^(O, Z a )4>{0, Z s ) dZ s 

Jr+xr™ Jr m « 

f%\t,Z s )d t <j>(t,Z s )dZ s dt. 



Now let us compute 

TV 



V / Vi- V x J N (t, Z N )<t>{t, Z s )lx N ev° N dZ N = ( div Xjv (V N f N {t, Z N ))<j>{t, Z s )lx N ev % dZ N 

using Green's formula. The boundary of V S N is made of configurations with at least one pair 
satisfying 1 < i < s and s + 1 < j < N, with \xi — Xj\ = e. 
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Let us define, for any couple £ [1, N] 2 , 

E^jl:=(z«eR M , \ Xi -Xj\=e 
(4.2.2) 1 

and V(M) e [l,a] x [s + 1, N}\ {i,j}, \x k - x t \ > e j . 

We notice that E^-(i, j) is a submanifold of {Zjy G R 2dAr , |xj — Xj\ = e], which is a smooth, codi- 
mension 1 manifold of R 2dN (locally isomorphic to the space x R d ( 2JV_1 )), and we denote by do % $ 
its surface measure, induced by the Lebesgue measure. Configurations with more than one collisional 
pair, i.e., and with 1 < i, i' < s, s + 1 < j,j' < N, with \xi — Xj\ = \xv — Xj>\ = e, 

and {i, j} ^ are subsets of submanifols of R dN of dimension at least two, and therefore con- 

tribute nothing to the boundary terms. 

X ' — x ■ 

Denoting := -y— — we therefore obtain by Green's formula: 



52 / Vi- V x J N {t, Z N )<f>(t, Z s )1 Xn £V% dZ N dt 

~[jR + xR 2dN 

= - V / f N (t, Z N ) Vi ■ V Xi (p(t, Z s )\ XNeV s dZ N dt 

i s N r 

+ ^12 52 / » hJ ■{vj-Vi)f N (t,Z N )<t>{t,Z a )da l gdt. 

i=1 j=s+1 Jn+x-z° N (i,j) 



By symmetry (1.1.1), this gives 



52 / v i -V x J N {t,Z N )<j)(t,Z s )l.x N ev%dZ N dt 

= - 52 / /jv(i, ^jv)vi • V Xi ^(t, Z s )l Xjv£l ^ dZ w dt 

+ A_(7V- S )V / ^ s+1 ■{v s+1 -v l )f N {t,Z N )(t ) {t 1 Z s )d ( j^ +1 dt ) 

v2 I=1 iR + xE5,(,, s +l) 



so finally by definition of f$ , we obtain 

N 

Vi ■ V Xi f N (t, Z N )(j>(t, Z s )lx N ev° dZ N dt 



(4.2.3) =-E/ ffl(t,Z s )v i -V Xi cf>(t,Zs)dZ s dt 

i=l JR+xR M » 

+ i__(iV- S )V / ^^+ 1 -(i; s+1 -« 4 )/ 7 v(t,Z Ar )^,Z s )da^ +1 dt. 

v2 i=1 -/R+xE^s+l) 



<R 2dJV 



4.3. CLUSTERS 
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Now let us consider the contribution of the potential in (4.2.1). We split the sum as follows: 

f F (^^A ■ VvMt, z N )<t>(t, z s )\ XNeV s dz N dt 
= - V / f f^i^i) . v Vi f N (t, z N )4>(t, z s )i XNeV s dz N dt 



N 



TLx N ev% dZ N dt . 

We notice that the second term in the right-hand side vanishes identically. It follows that 



-if 



F ( ) • V v Mt,Z N )<f>(t,Z a )l XireD . dZ N dt 



e 

so in the end we obtain 



Cl Xj > -v Vi 4>(t,z s )f^ } (t,z s )dz s dt 



f f^\t,Z s )(d t <f> + div Xs (V s <j>) + - F (— — —\ ■ S7 Vi <j>) (t, Z s ) dZ s dt 

JR+xR M » V £ ~~i V e / 7 

(4-2.4) = _f J^\0,Z s )<P(0,Z s )dZ s 

- E / " hS+1 ■ - "*) M*> z "W(t, z s ) dS N s+1 dt . 

V 1 i= i JR+x^ N (i,s+l) 

Remark 4-2.1. — Using the weak form of Liouville's equation, we see that multiple collisions (which 
occur as a boundary integral on a zero measure subset of dT> s N ) can be neglected. 



4.3. Clusters 

We want to analyze the second term on the right-hand side of (4.2.4). We notice that in the space- 
velocity integration the variables z s +2, ■ ■ ■ , zn are integrated over R d ( Ar_s_1 ) (with the restriction 
that they must be at a distance at least e from X s ) whereas z s+ \ must lie in the sphere centered 
at Xi and of radius e. It is therefore natural to try to express that contribution in terms of the 
marginal (Z s+ i). However as pointed out in (4.1.6), 

J il s+1) (Z s+1 )rfz s+1 ^./l s) (Z s ). 

The difference between those two terms is that on the one hand 

VX/v € T> , one has \xj — x s+ \ | > e for all j > s + 2 , 
which is not the case for Xm e T> S N , and on the other hand 

VX/v € V S N , one has \xj — x s+ \ \ > e for all j < s , 
a condition which does not appear in the definition of T> s ^ 1 . 

This leads to the following definition. 
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Definition 4-3.1 (e-closure). — Given a subset Xn — {x\,...,Xn} of R dN and an integer s 
in [l,N], the e-closure E(X s ,Xn) of X s in Xn is defined as the intersection of all subsets Y of Xn 
which contain X s and satisfy the separation condition 

(4.3.1) VyeF, VxeX N \Y, \x-y\>s. 
We denote \E(X S , Xn)\ the cardinal of E(X Si Xn). 

Now let us introduce the following notation, useful in situations where Xn belongs to Y, N (i,s + 1), 
defined in (4.2.2) 

Notation 4-3.2. — If X s+m = E(X S , X s+m ) and if for some integers jo < s < ko < s + m, there 
holds \xj — xk\ > e for all (j, k) £ [1, s] x [s + 1, s + m] \ {(jo, ko)}, then we say that E(X S , X s+m ) has 
a weak link at (jo,ko), and we denote X s+m = E(j ^ ) (X s , X s+m ). 

Moreover the following notion, following King [27], will turn out to be very useful. 

Definition 4-3.3 (Cluster). — A cluster of base X s = {xi,...,x s } and length m is any 
point {x s+ i, . . . , x s+m } in R dm such that E(X S , X s+m ) = X s+m . We denote A m (X s ) the set of 
all such clusters. 

The proof of the following lemma is completely elementary. 
Lemma 4-3-4- — The following equivalences hold, for m > 1 : 

(4.3.2) (e(X s , X n ) = X s+m ) (e(X s , X s+m ) = X s+m and X N E V s + m ) , 

/ T?(Y Y \ V \ / E {i,s+l)( X s,X s+m )=X s+m \ 

E(X S ,X N ) -X s+m \ x ^ V s+m\ 

\ \Xi - X s+1 \ = £ J 

as well as the implication, for m > 2, 

(4.3.4) [E(i tS+1 )(X s ,X s+m ) = X s+m ) ({x s+2 , ■ ■ ■ ,x s+m ) e A m _i(a; s+ i)J . 



4.4. Collision operators 

With the help of the notions introduced in Section 4.3, we now can reformulate the boundary integral 
in (4.2.4). 

Given 1 < s < N — 1 and Xn in T, N (i,s + 1), there holds |x s+ i — Xi\ = e, so that a; s+ i belongs 
to E(X S ,X N ), implying \E(X S ,X N )\ > s + l. We decompose T, s N (i,s + 1) into a disjoint union over 
the possible cardinals of the e-closure of X s in Xn ■ 

(4.4.1) E^(M + 1)= (J (x s N (i,s + l)f){Y N , \E(Y S ,Y N )\= s + m}), 

l<m<N-s 
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Xs 



Cluster of basis Xs+i 
And length m-1 



Xn\Xs 



Figure 4. Clusters with weak links 



implying 

i/< s+1 • (v s+1 - Vi ) f N {Z N )<j>{Z a ) d<r^ s+1 

= ^2 ^\E(x a ,x N )\=s+ m v l ' s+1 ■ {vs+Y-v.j) f N (Z N )4>(Z s )da l x S+1 . 

l<m<W- s "' S «(^ s + 1 ) 

By assumption of symmetry (1.1.1) for and <f>, if \E(X S , Xj^)\ = s + to, we can index the particles 
so that E(X a , Xw) = X s+rn : we obtain 



(4.4.2) 



TL\E(x s ,x N )\=s+ m v x,s+1 ■ {vs+i ~ Vi) f N (Z N )4>(Z s ) da 



i,s+l 
N 



c 



N-s-1 



^E(x 3 ,x N )=x s+ „y ■ {vs+i - Vi) fN{ZN)<P{Z s ) da% : 



i,s+l 



We use equivalence (4.3.3) from Lemma 4.3.4 and Fubini's theorem to write 

^E(x,,x N )=x s+ „y ,s+1 ■ {v s +i ~ v t )fN(Z N )<l>(Zs)da^ s+1 

= V2 [ ^ s+1 ■ (v s+1 ~ Vi)4>(Z a ) 

Js e (ii)xR J 

x ( / 1&, s+1) (x B ,x B+m )=x s+m fN +m Hz s+m )dZ( s+l s+m ) ) d<7i(x s+ i) , 

VJ R 2d (m -i) • y 

with do"i the surface measure on S E (Xi) := {a: € R d , |a; — £i| = e}. With (4.3.4), if m > 2, then the 
above integral over R, 2e2 ( m_1 ) appears as an integral over A m _i(x s+ i). We also remark that in the 
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case to — 1, we have a simple description of E^ is+1 - } (X S , X s+ i) = X s+ i : 

\Xi - X„+i I < £ 

>e for j e [1, s] \ {i} 



( \ ( \xi - x s+ i\ 

(4-4.3) (l, (i , a+1>M+0 ^ +1 ^o) <=► (^.^ 



This leads to the following definition of the collision term of order to > 1, for s + m < N : we define 

i) 



(4.4.4) 



C s . s+m /^ +m) (^ s ) := mC%_ s V / ^ S+M • (v s+1 - Vi 



x G <r s + 1 i>(/w +ro) )(^+ i ) <^i(a!«+i)dv 8 +i 



where for to = 1, by (4.4.3): 

(4A5) G<£ +1> ($ +1 >)(Z a+1 ):=( J] l ks+1 -„|> £ )/^ +1) (^ + i) , 



i<j< s 



and for m > 2 : 

-,(ro-l) / T(s+m) 



(4.4.6) y ~ (s+m) 

JA m _ 1 ( I « +1 )xR4".-i) 

The complex-looking indicator function T^E (i s+1) (x s .x s+m )=x s+m wm , in the estimates of the next 
chapters, be simply bounded from above by one. This will be the case for instance in an estimate 
showing that higher-order collision operators (4.4.6) are negligible in the thermodynamical limit; this 
estimate is (5.3.1) in Proposition 5.3.1. 

With (N - 3)0^21^ = mC N- s , we can now reformulate (4.2.4) into 

/ ffi(t, z s )(d t 4> + dw Xs (Vscf,) - 1 V v$ ( x i^ x l) . ( t , z s ) dz s dt 

Jr+xr 2 ^ v e t~Z x V £ J ' 

(4.4.7) 

r ~ N ~ s r ~ 

+ / /^(o, z a )4>(o, z„) dz s = Y] / </>(t, z s )c s , s+m f^ +m \t, z s ) dtdz s , 

JR M = m=1 JR+xR M ' 

so that fff appears as a (formal) weak solution to 

(4.4.8) djtf + V, 4 ft" - ' £ V$ • = £ C s , s+m / 



?"(s+m) 
' JV 



4.5. Mild solutions 

We now define the integral formulation of (4.4.8). Denote by & s (t) the s-particle Hamiltonian flow, 
and by H s the associated solution operator: 

(4.5.1) H.(i): /eC (n,;R)^/(*,(-i,.))eC fl (n,;R). 
The time-integrated form of equation (4.4.8) is 

(4.5.2) il s) (t,^) = H s (i)ft s) (0,Z s ) + V fu s (t-T)C s , s+m ^ +m \r,Z s )dT. 

m =i Jo 
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The total flow and total collision operators H and Cat are denned on finite sequences Gn = (g s )i<s<N 
as follows: 

Vs<N, (H(t)G N ) 6 :=H.{t)g s , 

N-s 

Vs < N — 1 , (C A rGjv) s := ^ C s ,s+m9s+m > (CatGat) n := . 



(4.5.3) 



We define mild solutions to the BBGKY hierarchy (4.5.2) to be solutions of 

(4.5.4) F JV (t) = H(t)F JV (0)+ / H^-TjCjvF^Tjdr, = (f^)i< s < N ■ 

Jo 

Remark 4-5.1. — At this stage, the use of weak formulations could seem a little bit suspicious since 
they are used essentially as a technical artifice to go from the Liouville equation (1.2.3) to the mild form 
of the BBGKY hierarchy (4.5.2). In particular, this allows to ignore pathological trajectories involving 
multiple collisions. Nevertheless, the existence of mild solutions to the BBGKY hierarchy (to be proved 
in the next two chapters) provides the existence of weak solutions to the BBGKY hierarchy, and in 
particular to the Liouville equation (which is nothing else than the last equation of the hierarchy). The 
classical uniqueness result for kinetic transport equations then implies that the object we consider, that 
is the family of truncated marginals, is uniquely determined (almost everywhere). 



Note that similarly we can define the total Boltzmann flow and collision operators S and C as follows: 

fV«>l, (S(t)G) s :=S s (t)g s , 
\V*>1, (C°G) s : C... ,.,/,.;. 
where S s denotes the free transport operator in s-particle space and C° s+1 is defined in (3.3.5). 



(4.5.5) 



CHAPTER 5 



CONTINUITY OF COLLISION OPERATORS 



In view of proving the existence of mild solutions to the BBGKY hierarchy (4.5.2), we need continuity 
estimates on the linear collision operators C s _ s+m defined in (4.4.4)-(4.4.5)-(4.4.6), and the total collision 
operator Cat defined in (4.5.3). 

We first note that, by definition, the operator C s . s+m involves only configurations with clusters of 
length m. Classical computations of statistical mechanics, presented in Section 5.1, show that the 
probability of finding such clusters is exponentially decreasing with m. 

It is then natural to introduce functional spaces encoding the decay with respect to energy and the 
growth with respect to the order of the marginal (see Section 5.2). In these appropriate functional 
spaces, we can establish uniform continuity estimates for the BBGKY (Section 5.3) as well as for the 
limiting Boltzmann collision operators (Section 5.4). 



5.1. Cluster estimates 

A point X s e R ds being given, we recall that A m (X s ) is the set of all clusters of base X s and length m 
(this notation is introduced in Definition 4.3.3 page 26). 

Lemma 5.1.1. — For any symmetric function ip on R Nd 7 any s <G [1,N — 1], any X s e R ds , the 
following identity holds: 

/ <p(XN)dX/ s+ i^ N \ = / lx N ev° <p(Xn) dX/ s+ i^ N \ 

J R (N-s)d J R d(N-s) 

(5 - L1) m f ff \ 

implying, for C > 0, 

(5.1.2) — f dX {s+hs+rn) < C m exp (Cn d (s + m)e d ) 
and 

C m+1 cxp(-CKd(m + l)e d ) f , . ... 

(5.1.3) I — / ^(2,m+i) <C(l-cxp(-C^ d )), 



m>l 

where Kd is the volume of the unit ball in R d 
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Proof. — The first identity (5.1.1) is obtained by a simple partitioning argument, which extends the 
splitting used to define C s . s+m in (4.4.4) in the previous chapter. We recall that, given any X s e R ds , 
the family 

|(ar s+ i, . . . , x N ) , \E(X S ,X N )\ = s + mj for < m < N - s , 
is a partition of H( N ~ s ^ d . Then we use the symmetry assumption, as we did in (4.4.2), to find 

/ <p(XN)dXf a+1 N ) = N Cn-s I ^E(x 3 .x N )=x s+m l p(XN)dX^ +1 _ N ' ) . 

It then suffices to use equivalence (4.3.2) from Lemma 4.3.4, noting that the set of all (x s+ i, . . . , x s+m ) 
in R md such that E(X s ,X s+m ) = X s+m coincides with A m (X s ). This proves (5.1.1). 

Estimates (5.1.2) and (5.1.3) come from the counterpart of (5.1.1) at the grand canonical level, i.e. 
when the activity is fixed, rather than the total number N of particles; Remark 5.2.3 expands on 
this terminology. 

For any bounded A C R d , the associated grand-canonical ensemble for n non-interacting particles is 
defined as the probability measure with density 

¥>n(X n ) := - P ^ - - Y\ ^x z eA- 

' l<i<n 

The s-point correlation function g s and the truncated .s-point correlation function g s are defined by 

°° n! f 
g s (X s ) :=^2 7Z ~r, / <fin(X n )dX (s+ltn) , 

°° n\ f 

g s {X s ) := 2_ \ 7 \T / lx„eBj^n(X„)(il( s+1 . n ) . 

n=s \ n — s >- JR(»-')i 

We compute 

f Vn(Xn)dX {s+hn) = C cxp ( - C|A|) (C|A| ]" S [] 



so that 



(C|A|)» 

Similarly, by definition of T) s n in (4.1.4), 



(5.1.4) 5s (X s ) = C s exp(-C|A|)^^J^ J] 1a(^)-C s Wa 

KKs 



/ ix n6 i?- n ix ie A^ (s+ i,„) = |An( fl c b £ (^))| = |An c s £ (x s )|, 

where we denote B e {X s ) := B £ (x i ) 1 with B £ (xi) := {y e R d , |y — < e}. This implies 



?.(x.)^^ ( -ciAi)E (aAn l B : ( y ir ' n- 



(n-s)! 



If B e (a;i) C A for all 1 < i < s, then |A| - |A n c 5 £ (A: s )| = \B £ (X S )\. We obtain 
(5.1.5) g s (X s ) = ( s cxp(-(\B £ (X s )\), UB £ (X s )cA. 



5.2. FUNCTIONAL SPACES 
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Besides, by (5.1.1), 
9s(X s ) = g s (X s ) 

°° n ~ s n\C m f ( f \ 

By Fubini, we get 

~ n \C™_ s f ( f \ 

n=sm=1 (« - S J ! ^A m (x s ) \7r("-3-™ )d / 

00 1 00 n! /" / /" \ 

oo If 

= E +1 (^I fc _ s( x/ fe(Xfe) ^ +1 ' fe) - 

We have proved that 

oo 1 . 

(5.1.6) ff a (X fl ) = 3s (X s ) + V — 9k(X k )dX (s+1 , k) . 

k=s+l i fc ~ s ) ! JA k ^(x s ) 

We now show how identities (5.1.4)-(5.1.5)-(5.1.6) imply the bounds (5.1.2)-(5.1.3). 

We first retain only the contribution of k = s + m in the right-hand side of (5.1.6). Given e > 
and X s e R ds , we choose A large enough so that B e {Y) C A for all Y G A m (X s ). This gives 

C > / C +m CXp ( - C|Be(X 5+m )|) dX (s + 1 . s+m) , 

and now \B e (X s+m )\ < KdS d (s + m) implies (5.1.2). 

We finally fix an integer K > 2 and choose s = 1 in (5.1.6). Given e > and x\ € R d , we choose A 
large enough so that B E (Y K ) C A for all Y K e Ajcfii). This gives 

K 

C-Cexp(-C|B e (a:i)|) >^ / exp ( - C|B e (*fe)l) dX (2>k) , 

fe=2 A fc _i(:zi) 

and bounding the volumes of balls from above, we find 

K ~ 1 /-fc+1 r 

C(l - cxp(-CK ( ie d )) > V ^T- exp ( - ( Kd (k + l)e d ) / dX {2 , k+l) . 
It then suffices to let K — » oo to find (5.1.3). This ends the proof of Lemma 5.1.1. □ 



5.2. Functional spaces 

To show the convergence of the series defining mild solutions (4.5.2) to the BBGKY hierarchy, we need 
to introduce some norms on the space of sequences (/^) s >i- Given e > 0, ft > 0, an integer s > 1, 
and a function g s : Q s — > R, we let 

(5.2.1) \9,\e,.,p~ SU P (\g.(Z 6 )\exp(pE e {Z a ))) 
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where for e > 0, the function E e is the s-particle Hamiltonian 

I 1 2 

(5.2.2) E e {Z s ) := ^ + E " *k) . with * e (x) := $ (*) , 

l<i<s l<i<fe<s 

and for e = 0, So is the free Hamiltonian: 

(5.2.3) E (Z S ):= ]T ^ • 

1<J<S 

Notation 5.2.1. — For e > and (3 > 0, we denote X E ^ s j the Banach space of continuous func- 
tions Q s — > R with finite \ ■ \ EtS ,/3 norm. 

By Assumption 1.2.1, for e > (and (3 > 0) there holds exp(/3F e (Z ;s )) — > oo as Z s approaches <9fi s . This 
implies for g s e X e<s ^ the existence of an extension by continuity: g s G C°(R 2ds ; R) such that g s = 
on dQ s , and g s = g on f2 s . 

For sequences of functions G = (g s ) s >i, with g s : f2 s — > R, we let for e > 0, /3 > 0, \i e R, 

IIGIIe.^.^ := sup ( | 5s | ejSi/3 exp(^s)) . 

s>l v y 

Notation 5.2.2. — For e > 0, j3 > 0, and /x e R, we denote X £i/ g ;jU i/ie Banach space of se- 
quences G = (g s ) s >i, with g s e ^ e , s „g and ||G|| £ii g i([1 < oo. 

The following inclusions hold: 

(5.2.4) if /?' < [3 and ^' < /j,, then X £i5;/3 / C X £)S ^ , X £ ^ iAJ , C X £j/3 ^ . 

Remark 5.2.3. — These norms are classical in statistical physics, where probability measures are 
called "ensembles" . 

At the canonical level, the ensemble e~^ Ee ^ Za ' > dZ s is a normalization of the Lebesgue measure, 
where f3 <~ T^ 1 (and T is the absolute temperature) specifies fluctuations of energy. The Boltzmann- 
Gibbs principle states that the average value of any quantity in the canonical ensemble is its equilibrium 
value at temperature T . 

The micro- canonical level consists in restrictions of the ensemble to energy surfaces. 

At the grand- canonical level the number of particles may vary, with variations indexed by chemical 
potential /ieR. 



5.3. Continuity estimates 

We now establish bounds, in the above defined functional spaces, for the collision operators defined 
in (4.4.4)-(4.4.6), and for the total collision operator Cat (4.5.3). In C s , s + m , the sum in i over [l,s] 
will imply a loss in /i, while the linear velocity factor will imply a loss in (3. The losses are materialized 
in (5.3.2) by inequalities /?' < /?, // < [i. 

Proposition 5.3.1. — Given (3 > and fi £ R, for to > 1 and 1 < s < N — to, the collision 
operators C SiS+m satisfy the bounds, for all Gn = (g s )i<s<N <= X £)J g )ft , 

(5.3.1) \C s>s+m g s+m {Z s )\ <e m - 1 C d e m ^{2^lf3) md l 2 [sr d/2 + ]T k|)e"^^ W ro U, s+roi/3 , 

l<i<s 
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for some Cd > depending only on d. 

If e < e^~ Kd (/3/27r) d//2 , then for all < f3' < f3 and // < ^t, the total collision operator Cjy satisfies the 
bound 



(5.3.2) \\C N G N \\ e ^^ < C d (l + /r d / 2 )(-^- + -J— ) ||G 

/or some C'd > depending only on d. 



N\\e,p,n , 



Considering the case to > 1 in (5.3.1), for which the upper bound is O(e), we see that higher-order 
interactions are negligible in the Boltzmann-Grad limit. 

Estimate (5.3.2), a continuity estimate with loss for the total collision operator Cat, is not directly 
used in the following. In the existence proof (Chapter 6), we use instead the pointwise bound (5.3.1). 

Proof. — We hrst consider the case m > 2. From the definition of G^™~^\ in (4.4.6), we obtain 

\G^ s +l\{g s +m){Zs+l)\ < \g s +m\e.s+m,l3 / exp ( — (3E e (Z s+m )) dZr s+2tS+m ) , 

where the norm | • | e . Sj( g is defined in (5.2.1), and the Hamiltonian E e is defined in (5.2.2). For the 
collision operator defined in (4.4.4), this implies the bound 

(5.3.3) \C S: s+rn,gs+m(Zs)\ < rnCN- s \gs+m\e,s+m,l3 * £] h,m ( Vs ) x Ji,m (X s ) , 

l<i<s 

where J, im is the velocity integral 

s+m 

/ / : 



/ (\v s +i\ + \vi\) exp ( - - \v J \ 2 )dV {s+hs+m) , 



and Ji.m is the spatial integral 

J t . m {X s ) := / exp(-/3 £ ®e( x j - x k)j da(x s+1 )dX {s+2 . s+m) 

JSe(x i )xA m - 1 (x. +1 ) l<j<k<s+m 

The velocity integral is a product of Gaussian integrals and can be exactly computed: 
(5.3.4) I hm (V s ) = (2n//3)( m - 1 W^(2n/^\v l \ + (2/f3) d ) exp ( - | £ \ Vj \ 2 ) . 

l<j<s 

For the spatial integral, there holds 



Ji, m {X s ) < exp ( - /3 £ ^e{xj - x k )j\S E (xi)\ x sup / dX ( 

V l<j<k<s X jA "-iW 

<exp(-/3 Y, *e{xj-x k j) xk^x ((m-^Ie^-^exp^Kd)), 



L (l,m-1) 



l<j'<fc<s 

where in the last bound we used identity (5.1.2) from Lemma 5.1.1 with s = 1 and ( = e~ d . This 
implies 

\C s , s+m9s+m (Z s )\ < C d e m - 1 ((N-s)s d ~ 1 ) m e m ^ {2k / P) md ' 2 (sp- d/2 + £ K|) 



Ki<s 



xe -PE e (Z e ) { I 

A e li/s+m |e,s+m,/3 • 



In the Boltzmann-Grad scaling Ne d 1 = 1, this gives (5.3.1). Above and in the following, Cd denotes 
a positive constant which depends only on d, and which may change from line to line. 
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In the case m = 1, by definition of in (4.4.5), there holds 

\ G< (!i + i)(9s+i)(Zs+i)\ < exp ( - PE s (Z a+1 ))\g a+1 \ £ta+1> p , 

and this implies 

\C SiS+1 g s+1 {Z s )\ < (TV- s)\g s+1 \ e , s+Uj exp (- (/3/2) ^ $ £ (^ - x fc )) x ^ /u(V s ) x k/" 1 , 

l<j'<fe<s l<i<s 

from which (5.3.1) is deduced as above. 

We turn to the proof of (5.3.2). From the pointwise inequality 

]T K| exp (-( 7 /2) ]T H 2 ) < S 1/2 (e 7 )- 1/2 , 7>0, 

l<i<s l<i<s 

we deduce for the above velocity integral Ii. m (V s ) the bound, for < j3' < f3, 

]T exp (0372) £ |^| 2 )^^V s )<C d (2^//3r d / 2 ( s r d/2 + S 1/2 (/3-/30- 1/2 )^ 

l<i<s l<J<s 

From the above bound in Ji, m (X s ), we deduce immediately, for < (3' < f3, 

max exp (> V & e (xj - x k ))j^ m (X s ) < n d (m - i)\ e m^ £ md-i 

Ki<s \ i — ' / 

l<j<fe<s 

With (5.3.3), these bounds yield, in the Boltzmann-Grad scaling, 

e ' E ' (Z ^' s \C^ s+m g s+m (Z s )\ < e m - 1 C d {^/P) md,2 e mK -e^ s { S p- d/2 + * 1/2 (/3 - PT 1 ' 2 ) 

X l^s+m I £,s+m,/3 • 

Summing over m, wc finally obtain, for Cat defined in (4.5.3), 

\\C N G N \\e,p>^ < C d \\G N \\ E .^ sup (( s p- d/2 + 5 1/2 (/3 - P')- 1 ' 2 )^-^ 8 ) 

l<s<N v ' 

x e- m ^- Kd) e m - 1 {2ir/l3) md/2 . 

l<m<N-s 

If e is small enough so that ee re<J_ '*(27r//?) d / 2 < 1, then the above series is convergent, and 

E -m(M-Kd) F rn-l ( o „ I Q\md/2 < e K " ^ (2n / f3) d/2 
l<m<N—s v ' ' 

Finally, 

sup f ( S f3- d/2 + s 1/2 (P - /3T 1/2 )e- (M -^ )s ) < e~\l + /3-" /2 )(m - m')" 1 + {P ~ PT\ 

\<s<N ^ > 

and this yields (5.3.2). Proposition 5.3.1 is proved. □ 

Remark 5.3.1. — We do not use the extra decay provided by the contribution of the potential in 
the exponential of the Hamiltonian. This is quite obvious in the bound for Ji, m {X s ) in the proof 
of Proposition 5.3.1, where we bound e ~^^ li a< k i'>+ m "M^-^) jjy e -0J2i< j<k <a "M^-x*^ Then, we 
might be tempted to replace E £ by the free Hamiltonian E in the definition of the functional spaces. 
The kinetic energy, however, is not a conserved quantity, so that in X . s _(s spaces the conservation of 
energy (6.1.5) does not hold. 
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5.4. Continuity estimate for the limiting collision operator 

Similarly to Chapter 2, we can define a limiting collision operator (see in particular (3.3.5) introduced 

in Chapter 3): 

(5.4.1) 

C° s . s+1 f (s+1) (Z s ) := W l,.(„. +1 - W4 )>o " ■ K+i - Vi ) 

x ^/( s + 1 )(f,a;i,wi, . . .,x h v*,. . .,x s ,v s ,Xi,v* +1 ) - f {s+1) (Z s , x u v s+1 )^dudv s+1 , 

where v* and v* +1 are obtained from Vi, v s+ i and i/ by the inverse scattering operator a^ 1 introduced 
in Chapter 3. The continuity estimate is as follows: 

Proposition 5.4-1- — Given (3 > 0, fi € R, the collision operator C® s+1 satisfies the following bound, 
for all g s+1 G X 0i8+li p : 

(5.4.2) \Cl a+l9s+1 (Z s )\<C d f3- d/2 (sf3- d/2 + £ |t>i|)e-^< z 'Wi|o,.+i,i8 , 

1<1<S 

/or some Cd > depending only on d. 
Proof. — There holds 

\Cs, s +i9s+i(Z s )\ < V / (\v s+1 \ + \v i \)(\g s+1 (v*,v* s+1 )\ + \g s+1 (v i ,v a+1 )\)di'dv 8+1 , 

omitting most of the arguments of g s+ \ in the integrand. By definition of | • | o,s,/3 norms and conservation 
of energy (3.1.4), there holds 

\g s+1 (v*,v* s+1 )\ + \g s+1 ( Vi ,v s+1 )\ < ( e -^o(z;) + e -/»^.>)| ffa+1 | 0i/3 

= 2e-^( z ')| ff8+1 | , a+ i lj8) 
where Z* is identical to Z s except for Vi and v s+ i changed to v* and This gives 

\Cl s+1 g s+1 (Z s )\<C d \g s+1 \ 0M e-P E °^ ]T J i;1 (V s ), 

l<i<s 

borrowing notation from the proof of Proposition 5.3.1, and we conclude with (5.3.4). □ 



CHAPTER 6 



LOCAL-IN-TIME WELL-POSEDNESS FOR THE BBGKY AND 
BOLTZMANN HIERARCHIES 



We state and prove an existence and uniqueness result for mild solutions to the BBGKY hierarchy, 
defined in (4.5.4), which we reproduce here: 

(6.0.1) F N (t) = U(t)F N (0) + [ U(t-T)C N F N ( T )dT, F N = (/^i^v , 

Jo 

as well as for the limit Boltzmann hierarchy 

(6.0.2) F(t) = S(t)F(0)+ f S(t-r)C°F(T)dT, F = (/^)i< s , 

Jo 

where the limiting collision operator C° as well as the free-particle flow S(t) are defined in (4.5.5). 



6.1. Functional spaces and statement of the result 

Existence and uniqueness for (6.0.1) will take place in spaces of X £;) g iA1 -valued functions of time (see 
Notation 5.2.2 page 34), where the indices f3 and \x themselves depend on time: in the sequel we 
choose for simplicity a linear dependence in time, though other, decreasing functions of time could be 
chosen just as well. Such a time dependence on the parameters of the function spaces is a situation 
which occurs whenever continuity estimates involve a loss, which is the case here since the continuity 
estimates on the collision operators lead to a deterioration in the parameters j3 and /i. We refer to 
Section 6.3 for some comments. 

Definition 6.1.1. — Given T > 7 a positive function (3 and a real valued function fi defined on [0,T] 
we denote X ej a. M the space of functions G : t € [0,T] i-> G(t) = (g s (t))i< s <G X e/ 3( t ) iAt ( t ), such that for 
all Z s e R 2ds , the map t e [0, T] g s (t,Z s ) is measurable, and 

(6-1.3) \\G\\e,e,»~ sup ||G(t)|| ei/S ( t ) >M (t)<oo. 

0<t<T 



We define similarly 



|||G|||o l/3 , M := sup \\G(t)\\ 0Mt)Mt) . 

0<t<T 



The existence result for the BBGKY hierarchy we shall prove is the following. 
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Theorem 2 (Uniform existence for the BBGKY hiearchy). — Let f3 > and /j, E R be 
given. There are T > and A > such that (3 := B — XT > 0, as well as e > and C > 0, such 
that for all < e < So, defining /j, := fiQ — XT, any family of initial marginals i*jv(0) — (/^' > (0)) seN , 
in X £i/ 3 0iAi0 gives rise to a unique solution F^ft) — (f]y\t))i< s <N in X £iJ g jA1 to the BBGKY hierar- 
chy (6.0.1) in the Boltzmann-Grad scaling Ne^ 1 = 1. It satisfies the following bound: 

l\F N \\ e ,e,„<c\\F N (p)\\e,eo,*>- 

This is a uniform existence result, in the sense that the existence time T does not depend on the 
number of particles N, which of course is crucial in the perspective of the limit N — > oo. Note that 
actually the only assumption made is on bounds on the initial family of marginals. 

For fixed e > 0, the uniqueness statement in Theorem 2 allows to define a maximal existence time T*(e). 
However we expect supT*(e) to be attained at e = 0, which precludes the definition of a maximal 

£>0 

existence time for the e-dependent family of hierarchies. We can however give a uniform bound from 
below for an existence time in Theorem 2: the following result is a corollary of the proof of Theorem 2, 
its proof is provided at the end of Section 6.2. 

Corollary 6.1.2. — Given (3 > 0, fi E H, for some constant Cd > 0, given 

(6.1.4) T:=C d e^(l + 2(3 d /2 )- 1 max Be^{(3 - B) d 

pe[o,/3 J 

then for all < e < £o, the solution to the BBGKY hierarchy with data i*V(0) = (//v'(0)) s£N * 
belonging to X e ^ 0iMo is defined on [0,T]. 

Remark 6.1.3. — For d <C 8n, there holds max Be~^(3o — P) d — + o(l)), hence an existence 
time of the order of f3^ 2 . 

A similar existence result as Theorem 2 can be obtained for the Boltzmann hiearchy. 

Theorem 3 (Existence for the Boltzmann hiearchy). — Let f3 > and € R be given. Then 
with the same notation as Theorem 2, any family of initial marginals F(0) = (/^(0))i< s *= ~^-o,^ ,no 
gives rise to a unique solution F(t) — (f^ (t))i< s in X ,/3, M to the Boltzmann hierarchy (6.0.2). It 
satisfies the following bound: 

\\\F\\\ o ^<C\\F(0)\\ OJjaM . 

The proof of Theorems 2 and 3 is typical of analytical-type results, such as the classical Cauchy- 
Kowalevskaya theorem. We follow here Ukai's approach [40], which turns out to be remarkably short 
and self-contained; the different approach of Nirenberg [32] and Nishida [33] would allow for direct 
use of the loss estimate (5.3.2). Let us give the main steps of the proof: we start by noting that the 
conservation of energy for the s-particle flow is reflected in identities 

1 5 ^ \ils(t)gs\e,sj3 = \g a \e,s,P and ||H(t)Gjv||e,/3^ = ||GW||e,/9 lA( , 

\Ss(i)go, s \o,s,p = \gs\o,s,p and \\S(t)G \\oj3^ = \\Go\\o,0,n , 

for all f3 > 0, /j, e R, g s e X Eta> p, g , s e X 0>St p, G N = (g s )i< s <N € X £i/8;jU , G = (go,s)i<s € Xo^^, 
and all t > 0. Next assume that there is a constant c < 1 such that for £ small enough (depending 
on c, Bq, fio, X and T) there holds the bound 

(6.1.6) V0<e<£ , / H{t-t')C N G N {t')dt' < c \\G N \\ e ^ . 

JO 
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This estimate is the object of Lemma 6.2.2 below. Under (6.1.6), the linear operator 



£ : G N e X £i/3iM ^|t4 Gjv(t) - jf H(t - t')C N G N (t') dt'^j 



ex. 



is is linear continuous from X £j| a M to itself with norm strictly smaller than one. In particular, it 
is invertible in the Banach algebra £(X £j/ g iM ). Next given Fn(0) <G X £j( g 0iAI0 , by conservation of en- 
ergy (6.1.5), inclusions (5.2.4) and decay of t i->- o — Xt and t i->- ^ — At, there holds 

(t^U(t)F N (0)) eX £ 

Hence, there exists a unique solution Fn e X £iJ g ijU to £Fjy = H(-)Fjv(0), an equation which is equiva- 
lent to (6.0.1). 

The reasoning is identical for Theorem 3, replacing (6.1.6) by 



(6.1.7) 



[ S(t-t')C°G(t')dt' < c II 0i/3 ,„. 



The next section is devoted to the proof of (6.1.6) and (6.1.7). 



6.2. Continuity estimates 



As explained in the previous paragraph, we need to prove (6.1.6), and its counterpart (6.1.7) for the 
Boltzmann operators. Let us first prove a continuity estimate based on Proposition 5.3.1. 

Lemma 6.2.1. — Under the assumptions of Theorem 2, there holds the bound, for < t < T, 



(6.2.8) 



f U s (t-t')C s , s+1 g s+1 (t')dt' <c(A),W),A,T)|||Gj V || e , / j, #t , 

JO e,8,f) -\t 

for all Gn = (g s +i)i<s<N € X £;J g ift , with c(/3o, Hoi A, T) computed explicitly in (6.2.14) below. 
Proof. — Let us define 

(6.2.9) /3$(t) := [3 - At and ^(t) := /J, -\t, 

so that (3 — Po(T) and n = H${T). By conservation of energy (6.1.5), there holds the bound 

/ U(t-t')C S)S+1 g s+1 (t')dt' < sup / e^^ E ^\C s . s+1 g s+1 (t',Z s )\dt' . 

Jo e,a,Po(t) z s en 2ds Jo 

Estimate (5.3.1) from Proposition 5.3.1 gives 
e^ E ^\C s , s+1 g s+1 {t',Z s )\ 



< 



C de ^(2 7 r/^(t')) d/2 |5 s+ i(t')l e , s+ i,^( t ')( s (/3o A (i'))- <i/2 + £ My'''-* 



)E e (Z s ) 



Ki<s 



By definition of norms || • || £ ,,3 jM and ||| • ||| £j| a, M we have 

l5. + i(^)| e , a+ i,^(to^ e " ( ' +lW(t ' ) ll G! ^(* , )ll e ,^(t'),^(*') 
< e -(^o(0||| GjV || £;/3 ^. 



(6.2.10) 



42 CHAPTER 6. LOCAL-IN-TIME WELL-POSEDNESS FOR THE BBGKY AND BOLTZMANN HIERARCHIES 



The above bounds yield, since /3q and ^o are nonincreasing, 



(6.2.11) 



-t)(s+E e (Z„)) 



f U{t-t')C S!S+l9s+1 {t')dt' 
Jo 

< \\G N \\ e .p^C d e^-^ T \2^/^{T)) d l 2 sup / C(t' ,t, Z s ) dt' , 

z s eR 2ds Jo 

where, for t' < t, 

(6.2.12) C(t',t,Z s ) := (,s(/3 A (O)^ /2 + E 

Since 
(6.2.13) 

there holds finally 

< c([3 , Ho, X,T)l\G N \\\ Et0ili . 



Ki<s 



sup f C{t',t,Z s )dt' < ^(1 



1 



£, S ,/3*(t) 



e^o(t) /" H(t-i')C s , s+l5s+1 (t / )rfi 
Jo 

where, with a possible change of the constant C^, 
(6.2.14) c(/3o,Mo,A,T) := C d e^° ^ , A * (l 



A(^(T))"/2 V ($(r))^ 



The result follows. 



□ 



In the next lemma, the definition (6.2.14) of c provides directly (6.1.6). 

Lemma 6.2.2. — Under the assumptions of Theorem 2, and for e small enough (depending 
on Po, Ho, A and T) there holds the bound, for <t <T, 

(6.2.15) V0<e<£ , f H{t - t')C N G N (t') dt' < 2c(/3 , Mo, A, T)\\\G N \\\ £ ^ , 

Jo e,/3 -At,Mo-At 

for all Gn — (g s )i<s<N € X £!| g :M , , where c is defined in (6.2.14). 



Proof. — We follow closely the proof of Lemma 6.2.1. The difference is that here we take into account 
higher-order collision operators C StS + m , with m > 2. 

Using notation (6.2.9), Estimate (5.3.1) from Proposition 5.3.1 gives 

gA) ( ) e( s ) |C S s _|_ m g s + m (t ,^ s )| 

< £ ™-V de ^(2 7 r//3 o A (i0) md/2 |5 s+ro (t0le, s+m ^(*o(«(^o(* / ))- <i/2 + E H)e A( *'- t)W) 

KKs 

Using also (6.2.10) with s + 1 replaced by s + m, we get 
/ H(t-i / )Cj V G JV (t / )<ft , || 

Jo IU,/3 A (*)^(t) 



(6.2.16) 



< IHGjvlHe^f E °™) SU P / C(t,t',Z s )dt'. 
y i<m<N-s ' Z.G&*>Jo 



where C m := C d e m ~ x e mi - R ^^ ( T » {2tt / (3 x {T)) md / 2 , and C is defined in (6.2.12) and satisfies (6.2.13). 
Under the assumption that 



(6.2.17) 



e e Kd -' x ^ T \2n/^{T)) d l 2 < 1/2, 
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we find 

(6.2.18) £ C m < 2C d e-^ T \^{T))- d / 2 . 

l<m<N-s 

The upper bounds in (6.2.13) and (6.2.18) are independent of s, and their product is equal 
to 2c(/?0) Mo ; A, T). Taking the supremum in s in (6.2.16) then yields (6.2.15). □ 

The proof of the corresponding result (6.1.7) for the Boltzmann hiearchy is identical to the first part 
of the proof of Proposition 8.1.1, since the estimates for C° s+1 and C s ,s+i arc essentially identical 
(compare estimate (5.3.1) from Proposition 5.3.1 with estimate (5.4.2) from Proposition 5.4.1). 

Proof of Corollary 6.1.2. — Given T, we are looking for A > and e > such that (6.2.17) holds 
and, say 

(6.2.19) C d {2 + {Bo - XT)- d/2 )e-^ +XT (B - XT)- d/2 = ^ ■ 

Indeed, if such a A exists, we can then define e = (l/3)e" Kti+AIO " AT (27r/(^o - AT))" d/2 , and (6.2.17) 
holds. With 8 = XT e (0,/? ), condition (6.2.19) becomes 

d 1 1 + 2(A,-^/ 2 
> C d e»°(l + 28* /2 r 1 pe-P(l3 - P) d , 
and (6.1.4) follows. □ 

In particular, given an existence time T for the BBGKY hierarchy, in the sense of Theorem 2, then T 
is an existence time for the Boltzmann hierarchy (6.0.2). 



6.3. Some remarks on the strategy of proof 



The key in the proof of (6.1.6) is not to apply Minkowski's integral inequality, which would indeed 
lead here to 



JQ 



H(t-t')C N G N (t')dt' 



< 



by (6.1.5), and then to a divergent integral of the type 



|C,/vG,/v(i') 



dt' . 



f n{t-t')c N G N (t')dt' 

Jo 



1 



+ 



1 



dt' 



The difference is that by Minkowski the upper bound appears as the time integral of a supremum in s, 
while in the proof of Lemma 6.2.2, and hence of (6.1.6), the upper bound is a supremum in s of a time 
integral. 

As pointed out in Section 6.1, other proofs of Theorems 2 and 3 can be devised, using tools inspired by 
the proof of the Cauchy-Kowalevskaya theorem: we recall for instance the approaches of [32] and [33], 
as well as [31] and [27]. 



CHAPTER 7 



ADMISSIBLE INITIAL DATA AND MAIN RESULT 



We state here our main result, describing convergence of mild solutions to the BBGKY hierarchy (6.0.1) 
to mild solutions of the Boltzmann hierarchy (6.0.2). This result implies in particular Theorem 1 stated 
in the Introduction. Existence and uniqueness results for both hierarchies were previously given in 
Chapter 5, as Theorem 2 page 40 and Theorem 3 page 40. 

The first part of this chapter is devoted to a precise description of Boltzmann initial data which are 
admissible, i.e., which give rise to solutions for which the convergence result holds. This involves 
discussing the notion of "quasi-independence" mentioned in the Introduction, via a conditioning of the 
initial data. Then we state Theorem 4 and sketch the main steps of its proof. 



7.1. Quasi-independence 

In this paragraph we discuss the notion of "quasi-independent" initial data. We first define admissible 
Boltzmann initial data, meaning data which can be attained from BBGKY initial data (meaning 
bounded families of truncated marginals) by a limiting procedure, and then show how to "condition" 
the initial BBGKY initial data so as to converge towards admissible Boltzmann initial data. Finally 
we characterize admissible Boltzmann initial data. 

7.1.1. Admissible Boltzmann data. — Let us define admissible Boltzmann initial data. 

Definition 7.1.1 (Admissible Boltzmann data). — Admissible Boltzmann data are defined as 
families F — {fQ S ^) s >i, with each f^ nonnegative, integrable and continuous over Ct s , such that 



JR 2d 

and which are limits of BBGKY initial data F Q _ N = (fQ%)i< s <N € X £i/ 3 0iM in the following sense: it 
is assumed that 



(7.1.1) 




(7.1.2) 



sup H-Fo.JvlkA^o < oo , for some ft>0,/i o eR, as Ne d 1 = 1 , 



AT>1 



(7.1.3) 



and 




U>„ (X N )f^ (Z N )dZ {s+hN) , 1 < s < N , 



46 



CHAPTER 7. ADMISSIBLE INITIAL DATA AND MAIN RESULT 



and that the following convergence holds: 

(7.1.4) fff^ — > f^ , for each s , as N — > oo with Ne^ 1 = 1 , locally uniformly in Cl s . 
In this section we shall prove the following result. 

Proposition 7.1.1. — The set of admissible Boltzmann data, in the sense of Definition 7.1.1, is the 
set of families of marginals F n as in (7.1.1) satisfying a uniform bound \\Fo\\o,p ,w> < 00 ■ 

7.1.2. Conditioning. — We first consider "chaotic" configurations, corresponding to tensorized 
initial measures, or initial densities which are products of one-particle distributions: 

(7.1.5) ff s {Z s )= H f ( Zi ), l<s<N, 

l<i<s 

where /o is nonnegative, normalized, and belongs to some Xq,i^ space (see Definition 5.2.1 page 34): 

(7.1.6) jo>0, / f (z)dz = l, / 6^o,iA for some /? > . 

JR 2d 

An important observation is that for (/® s )i< s <at defined by (7.1.5), with f satisfying (7.1.6), there 
holds in general sup ||(/® s )i< s <A r lle.^,,^ = +oo, for all (3 > 0, yu e R. Indeed, the correction in the 

Af>l 

Hamiltonian due to the potential $ e produces errors of size 0(1) in s-particle configuration subdomains 
such that \xi — Xj\ < e. These subdomains are not asymptotically small, even in the thermodynamical 
limit Ne 1 *- 1 = 1. 

This calls for cancelling out the contribution of the potential, by consideration of 

(7.1.7) / ®^(Z JV ):=exp(-A) ]T * e (xi- Xj j) J] > Z n^^n- 

l<i<j<N l<i<N 

With this definition, there holds the identity 

( 7 - L8 ) \fo,N \s,N,p = \fo\o,N,0 o , 

where the norms | • \e,N,p an d | • \a.N,p arc defined page 33. Indeed, using the notation for the 
Hamiltonian introduced in Section 5.2 page 33, there holds 

|/«fiSW = sup e ^ ZN ) e -p E e( x N ,o) lfo{Zi)l= sup Y[ e^ 2 ^\f (z t )\, 

z «e RMl i< s <N z «e RMN i< t <N 

and the last term in the right-hand side above is equal to | ./o I o^at /3 • 

The property of normalization is then preserved by introduction of the partition function 
(7.1.9) Z N := / ff$(Z N ) dZ N , l<s<N, 

and the definition of conditioned datum built on f as Z^f®^. This operation is called conditioning 
on energy surfaces, and is a classical tool in statistical mechanics (see [18, 29, 30] for instance). 

The partition function defined in (7.1.9) satisfies the next result, which will be useful in the following. 

Lemma 7.1.2. — Given f satisfying (7.1.6), there holds for 1 < s < N the bound 

1 < Z^Zn-s < (l - SKdlfolL^L 1 ) \ 

in the scaling Ne^ 1 = 1, where I/oIl^l 1 denotes the L°°(R^, L 1 (R^)) norm of f , and denotes 
the volume of the unit ball in R d . 
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Proof. — From the trivial lower bound 

exp(-f3 $ e (xi - xj)) > l| Xi _ Xj |> £ , 

we deduce the lower bound 

Z s+1 > y R2d(3+i) cxp ( - PoE e (X s , 0)) (J] l| Xl _ x . +1 |> e ) f® (s+1) (Z s+1 ) dZ {w) , 

with E £ (X S , 0) = ^ $s( x i — Xj), in accordance with the definition of the Hamiltonian (5.2.2). By 

l<i<j<s 

Fubini, we have 



f exp(-(3 E £ (X s ,0))( J] l lxi _ Xs+ll>£ )f® {s+1) (Z s+1 )dZ {w) 

y R M(. + D v i<4 1 

= / (TT a |xi-x. + i|>e)/o(^+l)^ S +l ) fo,N( Z s) dZ (l,s) ■ 



Since 

/ ( II a bi-a; s + il>e)/o( z s+l) rfz s + l > l/olii - KdS£ d |/oU= , 



we deduce from the above, by nonnegativity of the lower bound 

Z s+ i > Z s (\f \ L i - Kdse d \f \ L oa L i) , 

implying by induction 

JV-l 

Z N > Z N - S [j C 1 - 3£ d Hd\h\L^L^) > Z N - S (l - £Kd|/oU~Li) S ' 
j=N-s 

where we used s < N and the scaling Ne 4 ^ 1 = 1. That proves the lemma. □ 



7.1.3. Characterization of admissible Boltzmann initial data. — The aim of this paragaph 
is to prove Proposition 7.1.1. 

Let us start by proving the following statement, which provides examples of admissible Boltzmann 
initial data, in terms of tensor products. 

Proposition 7.1.2. — Given fo satisfying (7.1.6), define as in (7.1.7), and let fj^ be a con- 
ditioned datum built on fo ■ 

(7-1.10) &--=Z^f^. 

Then, families (f^%)i< s <N of truncated marginals of f^J , as defined in (7.1.3), satisfy (7.1.2) for 
any [1q such that e A ' |/o|o,/3 < 1- Moreover the data Fq = (/(f s ) s >i is admissible Boltzmann initial 
data associated with F a ,N — {fo%)i< s <N ■ 

Proof. — In a first step, we prove that untruncatcd marginals F 0<N := (fo%) s<N satisfy uniform 
bounds. In a second step, we prove that untruncated marginals converge uniformly in fl s towards 
We finally prove that truncated marginals converge as well. 
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First step. The trivial bound 



yields, using notation (7.1.7), 



Cxp(-A) E ^e(Xi-Xj))<l 

2<S,S + 1<_? 



&%{Z a ) < Z N \C N (Z S ) f exp ( - 13 V $ e ( Xi - Xj j) TT fo(zi) dZ (s+ 



1,JV) 

s+l<i<j'<Af ' s + l<i<N 

By symmetry, 

^AT— s : 



(7.1.11) / exp s+l,JV) — 2j\ 

jR 2d < N - s ) ' s +l<i<j<Ar s+l<i<N 

and this gives 



f ( s ) <-- <7 — 1 <7 

/0,JV - iV ^N-sJo^ 

< (l - £K d |/o|L~Li) S /^w , 

the second inequality by Lemma 7.1.2. 

By 2a; + ln(l - x) > for x e [0, 1/2], there holds 

(7.1.12) (1 - e^L/ok-LO^ < e 2s£K ^°l— 1 , if 2^1/oUoo^ < 1 , 
so that for N larger than some -/V (equivalently, for e small enough), 

spo|f( s ) I <■ p s(tto+2sKd\fo\ L co L i)\f®a I _ / „Mo+2eKd|/o| £ oo ii I f„ L „ V 

e l^o,iv| £il a ^ e Mo.jvU^o - ^ e l/0|0,/3 J ) 

the equality by (7.1.8). Given /z such that e Mo |/o|o,/3 < 1> f° r ^ larger than some iVi, which we may 
assume to be larger than N , there holds e fJ -° +2£Kd ^°^°° n |/ |o,/3 < 1. The above then implies 

sup H-Fiyvlk/Wo < sup sup fe^+^-^l^l/olo^^oo, 

AT>JVi N>Ni l<s<N v 7 

which of course implies the uniform bound sup ||fo,Af||e,/3 ,/j < °o- 

7V>1 

Second step. We compute for s < N : 

& = Z^ff; N l exp(- ]T /3b* c («i-^-)- E A><M*i-*i)) 

" /lt s+l<i<j<JV i<s<s+l<j 

x ] | /o(z») ^(s+l,AT) j 
s+l<i<N 

and deduce, by the symmetry property (7.1.11), 

(7.1.13) /q ^ = f® s N (^Z N - s - ^( s+li jv)) 
with the notation 



- / R2d(N _ s) i 1 - ex p (- E /*>*-(*< - 

■'«•'• i<s<s+l<j 

X exp ( - E Po®e(Xi-Xj)^ Yl h{ z i) dZ {s+l,N) 



s+l<i<j<N s+l<i<N 



so that Z( s+ i TV) is a function of X s . 
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From there, the difference — /gjy decomposes as a sum: 

(7.1.14) - = (l - Z N 1 Z N _^ff s N + -Z^^s+i.JV)/®^ ■ 

By Lemma 7.1.2, there holds 1 — Z jv 1 Zjv_ s — » as TV — > oo, for fixed s. Since /® s is uniformly bounded 
in Cl s , this implies that the first term in the right-hand side of (7.1.14) tends to as TV — > oo, uniformly 
in Q s . Besides, by < 1 — exp ^ — /3 ^ ®e( x i — x j)^j < E ^\xi-Xj\<e i we bound 

i<s,s+l<j l<i<s 

s+l<j<N 

Z\ s+hN) < J R2d(N _ s) ( E ^l^-x.-Ke) 

l^i^s- 7 "- 1 ' s+l<j<N 

X CXp ( - Po^eiXi-Xj^J Yl fo(Zi)dZ(s+l,N)- 

s+l<i<j<N s+l<i<N 

Given 1 < i < s, there holds by symmetry and Fubini, 



/ ( E A-lzi-XjKe) CX P ( - E Po®e(Xi ~ Xjfj ]J /o(^i) rf ^(s+l,JV) 

jR2d(jv- s ) s+ i<j<}f s+l<i<j<N s+l<i<N 

<{N-s) / l| x ._ Xs+1 | <£ /o(^+i)^ s+ i 
Jn 2d 

x / exp(- V /3o* e (xi - Xj)) ]~| /o(^i)d^( s +2,Jv) 

17 k ' s+2<j<j<.ZV s+2<i<W 

= (N - s) / l| Xi _ x | <e / (;Z s+ l)c£;2 s+ l X Zat_ s _i , 

so that 

(7.1.15) < s (^ - s)s d K d \f \ L oo L iZ N _ s _ 1 , 

where |/o|i°° z, 1 denotes the i°°(R x , i 1 (R^)) norm of / . By Lemma 7.1.2, we obtain 

Z N 1Z \s+l,N) < es^dl/olL^L 1 (l - eKdl/olL^L 1 ) , 

and the upper bound tends to as N — > oo, for fixed s. This implies convergence to 0, uniformly in Cl s , 
of the second term in the right-hand side of (7.1.14). 

We thus proved the uniform convergence 

fol " fo.N -> in Sl„. Since exp(/?^ £ (Z s )) -+ 1 locally 
uniformly in fl s (not uniformly in Q s ), the convergence f®%j — > /,f s holds locally uniformly in O s . We 
conclude that converges locally uniformly to tensor products /® s in f2 s . 

TTiini step. The bound (7.1.2) is a direct consequence of the corresponding bound for F ,n, proved 
in the first step, since < /q ' n < /q jy. The fact that the truncated marginals converge is due to the 
Lemma 7.1.3, stated and proved below. 

By the normalization condition in (7.1.6), the tensor products are marginals: 

f ff^) {ZsjZs+l)dZs+1 = f®s {Zs) f fo{zs + l)dZs+1 = f®s {Zs) 

This verifies (7.1.1), and concludes the proof. □ 
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Lemma 7.1.3. — Given Fq^n — (fo%)i<s<N satisfying (7.1.2) and (7.1.3) from Definition 7.1.1, 
with associated family F Q _ N = (f^ n )i< s <n of untruncated marginals: 

(7.1.16) ti s) N (Z s ) = f f^{Z N )dZ [s+1 . N) , l<s<N, Z s eCt s , jffl=ti N i}, 
there holds the convergence 

f^lt — f^ s lj — > , for fixed s > 1 , as N — > oo with Ne^ 1 = 1 , uniformly in fi s . 

Proof. — We apply identity (5.1.1) from Lemma 5.1.1 to f^, and obtain after integration in the 
velocity variables 

( ) ~( ) f + 

(7.1.17) f Q N (Z s ) - f Q N (Z s ) — > Cn-s / fo.N (Z s +m)dZ( s+ltS+m ) . 

Then, denoting Co = sup ||-Fo,M||e,,3o,Mo: a finite number by assumption, from 

M>1 

/^^ m) (Z s+m ) < cxp ( - n Q (s + m) - f3 E £ (Z s+m ))C 

<cxp(- Mo (s + ™)-(/V2) \ v i\ 2 ) C o, 

l<i<s 

we deduce, first by integrating the velocity gaussians and then by using the cluster bound (5.1.2) in 
Lemma 5.1.1 with ( = e~ d , the bound 

/ f^ m \z s+m )dZ {s+l , s+m) < (27r//3o) md/2 e-^ (s+m) Co / dX (s+l . s+m) 

JA m (X s )xTl d ™ JA m (X s ) 

< m\(2Tr/M md/2 e md e {Kd -' i « )< - s+m) Co . 

If 2ee Kd -M(27r//3 ) d/2 < 1, then 

N-s 

C%_ s m!{2n/p ) md/2 e md e (Kd -^ )(s+m '> < ^ (2ee K ^° (2ir/ fo) d/2 ) m — ► 

m— 1 m>l 

as e — > 0, implying — fffjf — > for fixed s, uniformly in fl s . □ 

Remark 7.1.4- — We can reproduce the above proof in the case of a time- dependent family of bounded 
marginals, i.e., Fn € Xj.^.^, with sup H-F/vlllfr,^ < °°: the notation of Definition 6.1.1. This 

N>1 

gives uniform convergence to zero, in time t G [0,T] and in space X s G tt s , of the difference between 
truncated and untruncated marginals: fffi — fffi — > 0. 

We now give the generalization of Proposition 7.1.2 that will be useful in the proof of Proposition 7.1.1. 
Let V = V{£l\) be the set of continuous densities of probability in Q,± : 

(7.1.18) V = {hGC°(Cl 1 ;K), h>0, [ h(z)dz = 1} . 

Let 7T be a probability measure on V, such that, for some /3 > and some ^ € R, 

(7.1.19) supprrc {h E T , \h\ ,p < e"" } . 
Next we define 

(7.1.20) 7r<"> := / h® s dit(h). 
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In the case that tt = 5f , then (7.1.20) reduces to the tensor product (7.1.5)-(7.1.6). We let 

h% s := exp ( - p ^ $ e ( Xi - Xjfjh?' , 2 < a < N , h® 1 = h e V , 

l<i<j<s 

generalizing (7.1.7), and 

(7.1.21) Z N := f h% N (Z N )dZ N , heV, 
generalizing (7.1.9). 

The following result is an obvious generalization of Lemma 7.1.2. 

Lemma 7.1.5. — Given ir satisfying (7.1.19) and h e suppn, the family of partition functions Z s 
defined in (7.1.21) satisfies for 1 < s < N the bound 

1 < Z~ X Z N _ S < (1 - eC d e-^p- 1/2 Y s , 

where Cd depends only on d. 

The next statement generalizes Proposition 7.1.2. Its proof is an immediate extension of the proof of 
Proposition 7.1.2 thanks to the dominated convergence theorem, using the obvious bound h% s < e~ s ^°. 

Proposition 7.1.3. — Given ir satisfying (7.1.19), let 

(7.1.22) tt^ V) := / Z N 1 h% N d7r(h) . 

Jv 

Then, families {^^)i< s <n of truncated marginals of as defined in (7.1.3), satisfy 

(7-1.23) SUp \\(^ ( N)l<s<N\\e,l3 B ,t,' < 1 

JV>1 

for any (i' < /io, with f3o and no from (7.1.19). The data (n^) s >i, with defined in (7.1.20), is 
admissible Boltzmann initial data associated with that family. 

By Proposition 7.1.2, tensor products (/ ( f" s ) s >i, with /o satisfying (7.1.6), are admissible Boltzmann 
data. More generally, by Proposition 7.1.3, the convex hull of the set of tensor products, in the sense 
of (7.1.19)-(7.1.20), is included in the set of admissible Boltzmann data. We finally show the converse: 
all admissible Boltzmann data belong to the convex hull of tensor products. 

We first remark that given a Boltzmann datum F , and an associated BBGKY datum F ,n, there 
holds 

(7-1.24) WFoh.fh.n <°°. 

with (3o and /zo as in (7.1.2). 

Indeed, let Co = sup H-Fb.ArlkAj.A'o < 00 ■ Given s and Z s e Q s , for e small enough, <& £ { x i — x j) = 0. 

N>1 

Besides, by (7.1.4) there holds the pointwise convergence f [ %{Z s ) -+ fk s) (Z s ). Hence taking the 
limit e in the left-hand side of the inequality e s ^ + ^ E ^ z ^\f^ ] N {Z s )\ < C , we find (7.1.24). 

The Hewitt-Savage theorem reveals the specific role played by tensor products: the set of families F = 
(/o^)s>i °f marginals (7.1.1) satisfying the uniform bound (7.1.24) is the convex hull of tensorized 
initial data, as described in the following statement. 
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Proposition 7.1.4- — Given Fo = (/q )«>i & family of marginals (7.1.1) satisfying the uniform 
bound (7.1.24) with constants /3q > and fi € R, there exists a probability measure ir over the set V 
of continuous densities of probability over Q.\, defined in (7.1.18), with 

(7.1.25) supp tt c {g e V, \g\ ,p < } , 
such that the following representation holds: 

(7.1.26) = / g^dn(g) , s>l. 

Jv 

Proof. — Given a family Fq satisfying (7.1.1)-(7.1.24), the existence of tt satisfying (7.1.26) is granted 
by the Hewitt-Savage theorem [23]. The goal is then to prove the inclusion (7.1.25). Assume by 
contradiction that, for some a > 0, 

(7.1.27) ir(A a ) = Ka >0, where A a := {g E T(R 2d ), \g\ ,i,p > + a} . 
We then have by (7.1.26) 

if > / f JS dAg), 

hence by f^ < e _SAio ||F || o,^ ,mo: we infer that ||-Fb||o,/9 ,A*o — + ae^ a ) s , which cannot hold for 
some a > and all s, since 1 + ae^° > 1. Hence (7.1.27) docs not hold, which proves the result. □ 

Proof of Proposition 7.1.1. — We already observed in (7.1.24) that admissible Boltzmann data are 
bounded families of marginals. Conversely given a bounded family of marginals F , by Proposi- 
tion 7.1.4 representation (7.1.26) holds. Then, by Proposition 7.1.3, F is an admissible Boltzmann 
datum. This proves Proposition 7.1.1. □ 

Combining Propositions 7.1.1 and 7.1.4, we see that all admissible Boltzmann data are built on tensor 
products, in the sense that given an admissible Boltzmann datum, representation (7.1.26) holds for 
some tt satisfying (7.1.25). 



7.2. Main result: Convergence of the BBGKY hierarchy to the Boltzmann hierarchy 
7.2.1. Statement of the result. — 

Our main result is a weak convergence result, in the sense of convergence of observables, or averages 
with respect to the momentum variables. Moreover, since the marginals are defined in Cl s , we must 
also eliminate, in the convergence, the diagonals in physical space. Let us give a precise definition of 
the convergence we shall be considering. 

Definition 7.2.1 (Convergence). — Given a sequence (/i^)k s <jv of functions h s N € C°(f2 s ;R), 
a sequence (h s ) s >i of functions h s G C (f2 s ;R), we say that (h s N ) converges on average and locally 
uniformly outside the diagonals to (h s ), and we denote 

(h s N )i< s <N — > (h s )i< s , 
when for any fixed s, any test function ip s e C^°(R ds ; R), there holds 

I v . (h' N - h s ) (X s ) := / tps{Vs)(h s N -h a )(Z a )dV a — >0, as N -> oo , 

locally uniformly in {X s G R ds , x,- L ^ Xj for i ^ j} ■ 
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With regard to spatial variables, this notion of convergence is similar to the convergence in the sense 
of Chacon. 

We remark that local uniform convergence in H s implies convergence in the sense of Definition 7.2.1: 

Lemma 7.2.2. — Given (/^ } ) i<s<w a bounded sequence in X £jJ g ;M with the notation of Defini- 
tion 6.1.1, if fffi — > /( s ) for fixed s, uniformly in t € [0, T] and locally uniformly in Cl s , then 
f ( N ] f (s \ uniformly in t e [0,T]. 

Proof. — Let K s be compact in {X s £ R ds , X{ ^ Xj for i ^ j} . There holds 

-/ (S) )(*-)| ^ IIV-ll^(R-) su p \{f { N ] - f {s) )(x s ,v s )\ . 

y s esupp</p 3 

The set K s x supp</? s is compact in O s . Hence the above upper bound converges to as TV — > oo, 
uniformly in [0,T] x if s . □ 

We can now state our main result: 

Theorem 4 (Convergence). — Given a potential that satisfies Assumption 1.2.1 stated page 3, 
given F admissible Boltzmann datum associated with a family (Fo,jv)at>i of BBGKY data, denot- 
ing Fpf the unique mild solution to the BBGKY hierarchy (4.5.2) with initial datum Fq.n, given by 
Theorem 2, and F the unique mild solution to the Boltzmann hierarchy (6.0.2) with initial datum Fq, 
given by Theorem 3, there holds the convergence 

F N ^F, 

uniformly on [0, T], for any common existence time T > 0. 

In particular, if the initial data F 0iN is asymptotically tensorized, meaning that F = (/q S ^) s >i 

s 

with flf\t,Z s ) = J^[/o(ij z i) then the first marginal fffl converges to the solution f of the Boltz- 
i=i 

mann equation (1.3.1) with initial data f - 

Solutions to the Boltzmann hierarchy issued from tensorized initial data are themselves tensorized. For 
such data, the Boltzmann hierarchy then reduces to the nonlinear Boltzmann equation (1.3.1), and 
Theorem 4 describes an asymptotic form of propagation of chaos, in the sense that an initial property 
of independence is propagated in time, in the thermodynamical limit. This corresponds to Theorem 1 
stated in the Introduction. 

The results in this chapter imply the following Corollary to Theorem 4. 

Corollary 7.2.3. Let F be a family of marginals (7.1.1) satisfying the uniform bound (7.1.24), 
and F be the solution to the Boltzmann hiearchy issued from F , as given in Theorem 3. There exists a 
family of solutions Fjv to the BBGKY hierarchy and Fn an associated family of untruncated marginals, 
such that 

F N ^ F and F N F , 
uniformly on [0,T], for any common existence time T. 
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Proof. — By Proposition 7.1.1, the family F is an admissible Boltzmann datum. Denoting -Fb.w 
an associated BBGKY datum, let T > be an existence time for the BBGKY hierarchy Fn with 
datum Fg jy, given by Theorem 2. By Theorem 4 the convergence I Vs ( fff — /^) — > holds uniformly 
in [0, T] and locally uniformly in O s . 

Then, by Lemma 7.1.3 and Remark 7.1.4, there holds fffi — fffi — > 0, for fixed s, uniformly in [0, T] x f2 s . 
By Lemma 7.2.2, this implies I Vs (f^ — /^) — > 0, uniformly in [0,T] and locally uniformly in fi s . 

We conclude that -A- f<- s \ uniformly in [0, T\. □ 

7.2.2. About the proof of Theorem 4: outline of chapters 8, 9 and 10. — 

The formal derivation presented in Chapter 2 (in the case of hard spheres, but which could easily be 
adapted to our case) fails because of a number of incorrect arguments: 

- Since mild solutions to the BBGKY hierarchy are defined by the Duhamel formula (4.5.2) where the 
solution itself occurs in the source term, we need some precise information on the convergence to 
take limits directly in (4.5.2). 

- The irreversibility inherent to the Boltzmann hierarchy appears in the limiting process as an arbitrary 
choice of the time direction (encoded in the distinction between pre-collisional and post-collisional 
particles), and more precisely as an arbitrary choice of the initial time, which is the only time for 
which one has a complete information on the family of marginals F .jy. This specificity of the initial 
time does not appear clearly in (4.5.2). 

- The heuristic argument which allows to neglect pathological trajectories requires to be quantified. 
These are 

- either trajectories for which the reduced dynamics with s-particles does not coincide with the 
free transport (H s ^ S s ), 

- or trajectories for which some of the (localized) interactions involve at least three particles (so 
that the scattering described in Chapter 3 does not apply). 

Indeed we have more or less to repeat the operation infinitely many times, since mild solutions are 
defined by a loop process. 

- Because of the conditioning by the energy E e (Zjy), the initial data are not so smooth. The operations 
such as infinitesimal translations on the arguments require therefore a careful treatment. 

To overcome the two first difficulties, the idea is to start from the iterated Duhamel formula, which 
allows to express any marginal f^\t, Z s ) in terms of the initial data Fo,n- By successive integrations 
in time, we have indeed the following representation of fffi: 

00 /■* /"tn-l 

H s+mi (ti — t 2 )C s+ 

mi .s+m2 ■ • ■ 

(7.2.1) „ =0 Mn Jo A) Jo 

...n s+mn (t n )f% +mn) (o) dt n ...dh 

where by convention f^\o) = for j > N, and the sum over M n :— (mi, . . . , m n ) is restricted to the 
range 1 < m, + i < TV — s — rrii with m := 0. 

Using a dominated convergence argument, we shall first reduce (in Chapter 8) to the study of a 
functional 

- which involves only the superdiagonal part of the collision operator (i.e. terms of the type Cjj + i), 

- defined as a finite sum of terms (independent of N), 

- where the energies of the particles are assumed to be bounded, 
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- and where the collision times are supposed to be well separated (namely \tj — tj + \\ > 5). 

The reason for the two last assumptions is essentially technical, and will appear more clearly in the 
next step. 

The heart of the proof, in Chapter 9, is then to prove the term by term convergence, dealing with 
pathological trajectories. Let us recall that each collision term is defined as an integral with respect 
to positions and velocities. The main idea consists then in proving that we cannot build pathological 
trajectories if we exclude at each step a small domain of integration. The explicit construction of this 
"bad set" lies on 

- a very simple geometrical lemma which ensures that two particles of size e will not collide in the 
future provided that their relative velocity does not belong to a small subset of R d (see Lemma 9.1.3), 

- scattering estimates which tell us how these properties of the transport are modified when a particle 
is deviated by a collision (see Lemma 9.1.4). 

This construction, which is the technical part of the proof, will be detailed in Chapter 9. The conclusion 
of the convergence proof is presented in Chapter 10. 



CHAPTER 8 



REDUCTIONS VIA DOMINATED CONVERGENCE 



The goal of this chapter is to use dominated convergence arguments to reduce the proof of Theorem 4 
to the term-by-term convergence of some functionals involving a finite (uniformly bounded) number of 
marginals (Sections 8.1 and 8.2). In order to further simplify the convergence analysis, we shall modify 
these functionals by eliminating some small domains of integration in phase space corresponding to 
pathological dynamics, namely large energies in Section 8.3 and clusters of collision times in Section 8.4. 

We consider therefore families of initial data: Boltzmann initial data Fq — (/o^) s eN such that 

|||-Pb|||o^o,w = sup sup (exp((3 E(Z s ) +^os)/ (s) ( 

sGN Z s 

and for each N, BBGKY initial data Fjv,o = (Jn\)i<s<n such that 



sup |||-Fjv,o||U,/3o,a*o = SU P SU P sup (exp(p E E (Z s ) + nos)ffi' (Z s )) < +oo . 

N N s<N Z s 

We then associate the respective unique mild solutions (constructed in Theorems 2 and 3 in Chapter 6) 
of the hierarchies 

/«(*) = S a (t)/ W + rS s (t-r)CL +1 / (s+m) (r)dr 
Jo 



and 

N— 



* S+m) (T)dT. 



(8.0.1) iS s) (t)=H s (t)/<; ) 0+ ^ / U s (t - r)C s , s+m fl 

In terms of the initial datum, they can be rewritten 

/ (s) (*) = e f r ■ ■ ■ r 1 s s(t tix.+is.+i^ - t 2 )c° s+hs+2 . . . 

„_ n Jo Jo Jo 



n=0 



...S s+n (t n )f^+ m ^(0)dt n ...dt 1 



and 



°° c-t rt\ t>t n — i 

In (*) = /> Z^, l I "'I H s (i — ti)C S;S+mi H s+TOl (t\ — t2)C s+mitS+m2 . . . 

n=Om u ...,m r Jo J ° J ° 



■ ■ ■ H s+mn (in)/^ +m "' ) (0) dt n . . . dt\ 



The observables we are interested in therefore involve infinite sums 

- because of multiple collisions in the BBGKY hierarchy, i.e. of collision terms C s , s + m where m is not 
uniformly bounded, 
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- and because there may be infinitely many collision times (n is unbounded). 



8.1. Reduction to first-order collision terms for the BBGKY hierarchy 



We will first prove that the estimates obtained in Chapters 5 and 6 enable us to reduce the study of 
the BBGKY hierarchy to the equation 

(8.1.2) g^(t,Z s ) = H s (t)f^(0,Z s )+ f H s (t - r)C s , s+1 ^ +1) (r, Z s ) dr , 1 < a < N - 1 . 

Jo 

Estimate (5.3.1) in Proposition 5.3.1 shows indeed that higher-order collisions are negligible in the 
Boltzmann-Grad limit. For the solution to the BBGKY hierarchy, this translates as follows. 

Proposition 8.1.1. — Let /3q > and /iq be given. Then with the same notation as Theorem 2, in the 



Boltzmann-Grad scaling Ne d 1 = 1, any family of initial marginals Fn(0) — ( Jn\Q)) 1 



<s<N 



€ X 



gives rise to a unique solution Gn G X^'J o ^ o of (8.1.2) and there holds the bound 



\\G n \\^t < C\\F N (0)\\ e ^ o . 

Besides, the solution Gn to the modified hierarchy (8.1.2) is asymptotically close to the solution Fn 
to the BBGKY hierarchy (6.0.1); 



(8.1.3) 

for some C > 0. 



Gn — Fn I 



< eC\\F N {0)\\e,(h,. 



In particular, given an existence time T for the BBGKY hierarchy, in the sense of Theorem 2, then T 
is an existence time for the modified hierarchy (8.1.2), in the sense of Proposition 8.1.1. 



Proof. — From Lemma 6.2.1, we deduce the existence and uniqueness result for (8.1.2), and the bound 
for Gn, in the same way that Lemma 6.2.2 implies Theorem 2. Notice that an existence time for the 
BBGKY hierarchy is an existence time for the modified hierarchy, since the bound (6.2.8) is better 
than (6.1.6). 



We turn to the proof of (8.1.3). There holds 



\\G 



• ~ FnW^.t 



< 



H s (t-t')C s , s+ i(^ +1) -./l S+1) )(t') 



dt' 



Ks<N 



+ 



f{a s (t-t') ]T c s , +m f^ +m) (t')) i 



dt' 



<s<N 



With (6.2.8), this implies 

\\Gn — Fn\\ x *,t < cq 



/*(H.(t-f) £ c s , s+m ti +m \t')) 

J ° 2<m<N-s 



dt' 



<s<N 



with Co := (l — c(/?o, fJ-o, A, 7 1 )) , which is indeed strictly positive by assumption. We conclude as in 
the proof of Lemma 6.2.2. □ 



From now on we therefore shall concentrate on equation (8.1.2) instead of (8.0.1). 



8.2. REDUCTION TO A FINITE NUMBER OF COLLISION TIMES 
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8.2. Reduction to a finite number of collision times 



By successive integrations in time of (8.1.2), we get a representation of gff in terms of the initial 
datum q , for all n such that s + n < N : 

00 ft r t\ rt„ 



(8.2.1) 



n=0 



where by convention f^ = for j > N 



o Jo Jo 



H s (i — ti)C s , s+ iH s+ i(ti — t 2 )C s+ i^ s+ 2 ■ ■ ■ 



■ ■■n s+n (t n )f^ +n) (0)dt n ...dt 1 



Similarly for mild solutions to the Boltzmann hierarchy, we have 

/ (s) w=£ f r ■■ r i s S (i-ti)c s . s+1 s s+1 (i 1 -t 2 )c s+1 , s+2 ... 

n=0 Jo Jo Jo 

...S s+ „(t„)/< s+n) (0) dt n ...dh. 



(8.2.2) 



Due to the uniform bounds derived in Chapter 6, the dominated convergence theorem implies that it 
is enough to consider finite sums 

9N n) (t) = E I I 1 ■ ■ ■ I " 1 H ^ - *i)C 8 , a +iH a+ i(ti - t 2 )C s+M+2 . . . 
k=o 



(8.2.3) 



" r t r ti 
k~0 J ° J ° 



...H s+k (t k )f^ +k) (0) dt k ...dt x 



S s (i — ti)C s , s+ iS s+ i(ti — t2)C s +i, s +2 ■ ■ ■ 
■■■Ss + k(t k )ti s+k) (0) dt k ...dh 



and therefore to study the term- by-term convergence (for any fixed k), as expressed by the following 
statement. 



Proposition 8.2.1. — The following estimates 



< C [ ~- ) \\F Nfi \\ e ^ e s ^ 



■ <"( = ) U^lk/W^ 



f (s) {t) _ f{ s,n) {t) 

hold uniformly in Z S ,N and t <T, in the Boltzmann- Grad scaling Ne^ 1 = 1. 



Proof. — Using the continuity estimate (6.1.6) together with the condition (6.2.19) on A, we get 



(8.2.4) 

We then deduce that 



f U(t-t')C N G N {t')dt' 
Jo 



S,0O — At,|Uo — Xt O e,/3 ,(i 



O s(n -Xt) 



E 



5.2.5) 



k=n+l 



'1 f^k-l 



o Jo Jo 



H s (t — ti)C s , s +iH s +i(ti — i2)Cs+i,«+2 • • • 



...H 8+fc (t fc )^r' +fc) (0)dtfc...dti 



< C . 

l~ ~ V 3 



|Gjv|| x a,t 
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Combining this estimate together with the uniform bound on ||Gjv|| x A - T leads to the first statement 
in Proposition 8.2.1. 

The second statement is established exactly in an analogous way, using estimate (5.4.2) together with 
the uniform bound obtained in Theorem 3. □ 

From now on we therefore consider the approximate observables : 

(8.2.6) I s , n (t)(X s ) := I ' y s (V s )t°' n \t,Z s )dV s , and I° s Jt)(X s ) := I \p s (V s )f^ n \t, Z s )dV s . 



8.3. Energy truncation 

Recall the definitions of the Hamiltonians given in (5.2.2) and (5.2.3): 

I |2 I 

E £ (Z S ):= * e (*i - *fc) with $ e (a;) := * (^) , and E (Z S ) := ]T ^ 



l<i<s l<i<k<s 

We introduce a parameter R > and define 



Ki<s 



= £ / / 1 • • • / ^ H a (t - t 1 )C a , a+ iH a+1 
k=o 



(ti - * 2 )C. 



S+1,8+2 ■ 



■ ■'H- s+k (tk)TL\E E (z B+k )\<R 2 fN +k \o) dt k ...dh, 



1.3.1) 



— /\, / / ••• / S s (i — ti)C s , s+ iS s+ i(ti — t 2 )C s+ i !;s+ 2 . . . 
k=Q Jo Jo Jo 

■ ■ ■ S s+ fc(ife)l|B (z s+fc )l<« 2 /o" +fe) (0) ■ • ■ d*i 

and the corresponding observables 

(8.3.2) := J Vs {V s )g { °${t,Z s )dV s , and J?,„(t)(X 8 ) := J ^ s (V s )f R s ' n) dV s . 

Using the bounds derived in Chapters 5 and 6 we find easily that I s , n {t) ~ I R n (t) anc l ^s.n W — ^s.'jf (*) 
can be made arbitrarily small when R is large. More precisely the following result holds. 

Proposition 8.3.1. — There are some nonnegative constants C,C depending only on (s,n,t) such 
that 



\\Is,n(t, J, M) — I^ n (t, J, M)\\ L oc( R d8) < CIMIioo^e- R \\\F N ,o\\\e,p ,»o , 



and 



|/ s °„(i, J,M) -I° s £(t, J.AOIU^r-.) < Cyhoo^e-c'^WlFoWo,^ 



Ho ■ 



Proof. — Let (3' < (3 be such that (3' — Xt > 0. The estimates in Chapter 6 show that 
tR (+\<v \\ s n ii,„ii . . . o„„„„~ n. _ ._ ... _7(s+k) 



\I B>n (t)(X e ) -I s n (t)(X s )\ < C„||(^|| L oo (R£!s) SUpSUp l\E e (Z e+k )\>R*fN,0 

r<t k<n 



s,s-\-k,j3' — cr 



<Cy s \\ LoaiRd s ) e^-M R2 11^(0)11^^ 
which proves the result, with C ~ /?o — ct. The argument is identical for I® n (t) — I^{t). 



□ 



8.5. REFORMULATION IN TERMS OF PSEUDO-TRAJECTORIES 
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Remark 8.3.1. — It is useful to notice that the collision operators preserve the bound on high ener- 
gies, in the sense that 

C s ,s+1^E E (Z S+1 )<R 2 = ^-E E (Z S )<R 2 Cs,s+l^-E e {Z s+1 )<R? 
C s ° s+ ll£(Z J+1 )<ff = ^E(Z e )<R 2 Cl s+1 t E (z s+1 )<R 2 ■ 



8.4. Time separation 

We choose another small parameter <5 > and further restrict the study to the case when ti — t i+1 > S. 
That is, we define 

T n (t) := |t„ = (ti, . . . ,t n ) I U < ti-i with t n+ i = and t = ij , 
T n ,s(t) :={T„eT /U-t i+1 >6}, 



and 



(8.4.1) 



I* n 5 (t)(X s ) := / tp B {V B ) [ H a (t - ti)C a , a+ iH a+ i(ii - t 2 )C s+w 
J JT n , s (t) 

■ ■ ■ Cs+n-l,s+n,H- s+n (t n — t n+ i)t\ Ec ( Zs+n )\<R 2 InJ)^ ' dT n dV s , 

I° s ;«' S (t)(X s ) := / Vs {V s ) f S s (t- ti)C? ii+1 S, +1 (ti - t 2 )C° s+w 
J JT n ,s(t) 

■ ■ ■^ a s +n-l,s+n^'s+n{tn ~ *n+l ) 1| E (Z s+n ) | <R 2 fo^^ dT n dV s . 



Again applying the continuity bounds for the transport and collision operators, the error on the 
functions I* n (t) - I*£(t) and (t) - can be estimated as follows. 

Proposition 8.4-1- — There is a constant C depending only on (s,n,t), such that 

II -^n(*) _ l£n (*)IU°°(R<»«) < C^I|y||z,«=(R.«*«) ID -PjV.O |||er,/3 ,/*o 

and 

Wlslnit) I° s ^ 5 (t)\\ L oo {n ^ < CJ|MU» (R ,. ) ||F ||o,/»b,Mo • 



8.5. Reformulation in terms of pseudo-trajectories 

In the integrand of the collision operators C s , s +i defined in (4.4.4), we now distinguish between pre- 
and post-collisional configurations, as we decompose 

C s ,s+1 = Ct,s + i ~ C ss+1 

where 

(8-5.1) C± s+1 / (s+1) = £ C± : ™J^ 

rn — 1 
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the index m referring to the index of the interaction particle among the s "fixed" particles, with the 
notation 

(C s ± ;+i/ (s+1) )(^) := (N- S )e 2 f (u-(v s+1 -v m )) ± J^ +1 \Z s ,x m + eu,v s+1 ) 

x n a 

\xj-x B+1 \>s dvdv s+ \ , 

i<j< s 

the index + corresponding to post-collisional configurations and the index — to pre-collisional config- 
urations, according to terminology set out in Chapter 3. 

In the same way we have to decompose the Boltzmann collision operators (5.4.1) into 

r o _ r o,+ _ r o- 

where the index + corresponding to post-collisional configurations and the index — to pre-collisional 
configurations. By definition of the collision cross-section in Chapter 3, we have 



(C° s \+Tf {s+1) )(Zs) := / KK+i - v m \,w)f(' +1 \Z e ,x m ,v B + 1 )dudv a+ i 

Js d 1 - 1 xR d 

= / ((v s+1 - v m ) ■ v) + f { - s+1) (Z s ,x m ,v s+ i)dvdv s+ i and 

(C°' s + +i7 (s+1) ) (Z s ) := / - «m|, w)/< s+1 >(*i, • • • , x m , v* m , . . . , z s ,x m , v* s+1 ) duidv s+1 

•/S? _1 xR d 

= / ((v s +i - v m ) ■ v) + f {s+1 \zi 1 . . . ,x m ,v* m , . . . , z s ,x m ,v* s+1 )dudv s+ i . 



lsf- 1 xR d 

f variable 
for the BBGKY collision operators. 



Performing the change of variables v ^ —v in the integral defining C° s s _^ , we get similar formulas as 



The BBGKY and Boltzmann observables we are interested in (see the notation of Definition 7.2.1) can 
therefore be decomposed as 



I%(t,X,) (f[ji)l^(t,J,M)(X s ) and 

n=0 J.M i=l 



n=0 J.M 

where the elementary junctionals /j^f(t, J, M) are defined by 



I^ s (t,J,M)(X s ) := / <p a (V a ) f H a (t-t 1 )C%™lH a+1 {t 1 -t 2 )C*ffi +2 
J JT n , s (t) 

■ ■ ■ H s+ „(t„ — t n -\-i)1\E e (z e+n )\<R 2 fNfi n ^dT n dV s 



5-5-3) , , 

I° s: y(t, J, M)(X S ) := / Vs (V s ) / S s (t - OC s °;fVr S s+ i(ti - t 2 )c°4i : ^ 2 



• • • S s+ „(i„ — t n+ i)\\ Eo (z s+n )\<R? fo S+n ' ] dT n dV s 

with 

J := (ji, . . . , j n ) e {+, -}™ and M := (mi, . . . , m n ) with € {1, . . . , s + i - 1} . 
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Each one of the previous functionals I^(t, J, M) and I^' 5 {t, J, M) defined in (8.4.1) can be viewed as 
the observable associated with some dynamics, which of course is not the actual dynamics in physical 
space since 

- the total number of particles is not conserved; 

- the distribution does even not remain nonnegative because of the sign of loss collision operators. 
This explains the terminology of "pseudo-trajectories" we choose to describe the process. 

In this formulation, the characteristics associated with the operators H s+i (ti — t i+ i) and S s+i (ti — t i+ i) 
are followed backwards in time between two consecutive times U+i and U, and collision terms (associated 
with C s+ i tS+ i + i and Cg +i s+i+1 ) are seen as source terms, in which, in the words of Lanford [31], 
"additional particles" are "adjoined" to the marginal. 

The main heuristic idea is that for the BBGKY hierarchy, in the time interval [ij+i,tj] between two 
collisions C s _|_j_i jS+ j and C s +j iS +i+i, the particles should not interact in general so trajectories should 
correspond to the free flow S s +i. On the other hand at a collision time U, the two particles in interaction 
may be subject to the scattering operator and thus their velocities are liable to change. This is depicted 
in Figure 5. 



i.(Wi) 

Zj(t„ + l) 



Zs+2(£ n +l) 
X, + J.tn+l) 



+1 = 

































x 2 












. Il(tl) 


Il(tj) 

















Figure 5. Pseudo-trajectories 

At this stage however, we still cannot study directly the convergence of I^(t, J, M) — I^^' S (t, J, M) 
since the transport operators H fc do not coincide everywhere with the free transport operators Sfc, 
which means - in terms of pseudo-trajectories - that there are recollisions. Note that, because the 
interaction potential is compactly supported, recollisions happen only for characteristics such that 
there exist i,j <E [1, k] with i ^ j, and r > such that 

\(xi - tv 1 ) - (xj - TVj)\ < e . 

We shall thus prove that these recollisions arise only for a few pathological pseudo-trajectories, which 
can be eliminated by additional truncations of the domains of integration. This is the goal of the next 
chapter. 



CHAPTER 9 
ELIMINATION OF RECOLLISIONS 



We have seen in the previous chapter that the convergence of observables (stated in Theorem 4 in 
Chapter 7) reduces to the convergence to zero of the functional lf-£ — I^^' 6 'i defined in (8.4.1), 
corresponding to dynamics 

- involving only a finite number s + nof particles, 

- with bounded energies (at most R 2 ), 

- such that the n additional particles are adjoined through binary collisions, 

- at times separated at least by S. 

Let us denote, for any constant c > 0, by Gk(c) the set of "good configurations" of k particles, separated 
by at least c through backwards transport: that is the set of (Xk, Vk) € R dfe x B R such that the image 
of (Xk, Vk) by the backward free transport satisfies the separation condition 

Vr > 0, Vz 7^ j, \xi - Xj 

in particular it is never collisional. We recall that B k R : 
we write Br := B R . 

Our goal in the present chapter is to slightly modify (in a uniform way) the functionals 1^ and Ig'^' 5 
in order for the corresponding BBGKY pseudo-trajectories to be decomposed as a succession of free 
transport and binary collisions, without any recollision. Paragraph 9.1 is devoted to the statement and 
the proof of a geometrical proposition showing how to eliminate bad sets in phase space so that for any 
particle outside such bad sets, adjoined to a good configuration, the resulting configuration is again a 
good configuration. This is applied to the Boltzmann and BBGKY hierarchies in Paragraph 9.2. 



- r(vi - Vj)\ > c, 

= {v k eR dk /\V k \<R] and in the following 



9.1. Elimination of bad sets in phase space leading to recollisions 

9.1.1. Statement of the result. — In this section we momentarily forget the BBGKY and Boltz- 
mann hierarchies, and focus on the study of pseudo-trajectories. More precisely our aim is to show 
that "good configurations" are stable by adjunction of a collisional particle provided that the deflection 
angle and the velocity of the additional particle do not belong to a small pathological set. Furthermore 
the set to be excluded can be chosen in a uniform way with respect to the initial positions of the 
particles in a small neighborhood of any fixed "good configuration" . 
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Notation 9.1.1. — In all the sequel, given two parameters r\\ and r\i, we shall say that 

m < m m < Cf] 2 

for some large constant C which does not depend on any parameter. 

Proposition 9.1.1. — Let a,e ,n <C 1 be such that 

(9.1.1) a<£ <?7(5. 

Given Z k G G k {e n ), there is a subset B k (Z k ) ofSf -1 x B R of small measure 



(9.1.2) \B k (Z k )\ <C(i?)r 7 d - 1 + C(i?,$,r/)( 



a_ + £o 

£ S 



(9.1.3) Vr>0, 



and such that good configurations close to Z k are stable by adjunction of a collisional particle close 
to x k and not belonging to B k (Z k ), in the following sense. 

Consider [y, v) G (S^ -1 x Br) \ B k (Z k ) and let Zk be a configuration of k particles such that Vk = V k 
and \Xk — Xk\ < a. 

• If v ■ (v — Vk) < then for all e > sufficiently small, 

Vi ^ j E [1, k] , \(xi- TVi) - (xj - TVj)\ > e , 

Vj G [1, k] , \(x k +ev- tv) - (xj - TVj)\ > e . 

Moreover after the time S, the k + 1 particles are in a good configuration: 
(9.1.4) {X k -5Vk,V k ,x k + eis-5v,v) eG k +i(£o/2). 

• If v ■ (v — Vk) > then define for j £ [1, k — 1] 

(4* ,z £ *) :=a~ 1 (z k ,(x k +su,v)) and zf := (xj - t s Vj, vj) , 

where a e is the scattering operator as in Definition 3.2.1 and where t £ denotes the scattering time 
between z k and (x k + ev, v). Then for all e > sufficiently small, 

' Vz ± j e [1, k] , \{xf - rvD - {xf - rvf )\ > e , 
Vj g [1, fc] , |(ar £ * - r«") - (xf* - rvf)\ > e . 

Moreover after the time 5, the k + 1 particles are in a good configuration: 

(9.1.6) (xr-(5-t £ )V k s *,Vt,x s *-(6-t s )v s *,v**) e g k+1 (e /2) . 

The proof of the proposition may be found in Section 9.1.3. It relies on two elementary geometrical 
lemmas, stated and proved in the next section. The first one describes the bad trajectories associated 
to the (free) transport. The second one explains how they are modified by the scattering. 

Remark 9.1.2. — For the sake of simplicity, we have assumed that the additional particle collides 
with the particle k. Of course, a simple symmetry argument shows that an analogous statement holds 
if the particle k + 1 is added close to any of the particles in Z k . 

The proof of Proposition 9.1.1 shows that if Z k = Z k then the factor e a /2 in (9.1.4) and (9.1.6) may be 
replaced by £q. The loss comes from the fact that the set to be excluded has to be chosen in a uniform 
way with respect to the initial positions of the particles in a small neighborhood of X k . 



(9.1.5) Vt>0 : 



9.1. ELIMINATION OF BAD SETS IN PHASE SPACE LEADING TO RECOLLISIONS 
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9.1.2. Two geometrical lemmas. — In this section we state and prove two geometrical lemmas 
which though elementary, are the key to the proof of Proposition 9.1.1. Here and in the sequel, we 
denote by K(y, p) the cylinder of axis y £ R d and radius p > and by B p (y) the ball centered at y of 
radius p. 

Lemma 9.1.3. — Given a > satisfying e<C«C sq, consider x\, x 2 in R d such that \x~i — x 2 \ > £o> 
and v\ € Br. Then for any x\ € B a {xi), any x 2 <E Ba(x 2 ) and any v 2 € Br, the following results 
hold. 

• If vi — v 2 $l K(xi — x 2 , 6Ra/e ), then 

Vr>0, |(cci - v-tr) - (x 2 - v 2 t)\ > e; 

• If vi - v 2 K(xi - x 2 , 6Re /S) 

Vr > S , \{xi- Vit) - (x 2 - v 2 t)\ > e . 

Proof. — • Assume that there exists such that 

|(xi - «it) - (x 2 - v 2 t)\ < e. 

Then, by the triangular inequality and provided that e is sufficiently small, 

\{x\ - x 2 ) - t*(vi - v 2 )\ < e + 2a < 3a . 

This means that (vi — v 2 ) belongs to the cone of vertex based on the ball centered at x\ — x 2 and of 
radius 3a, which is a cone of solid angle (3a/|a;i — a^Q^ 1 (since a <C £o)- 

The intersection of this cone and of the sphere of radius 2R is obviously embedded in the cylinder of 
axis X\ — x 2 and radius 6Ra/e , which proves the first result. 

• Similarly assume that there exists r* > S such that 

K^i - vit) - (x 2 - v 2 t)\ < e . 
Then, by the triangular inequality again, 

\(xx - x 2 ) - t*(vi - v 2 )\ < e + 2a < 3e • 
In particular, for any unit vector n orthogonal to xq — x, 

r*\n ■ (vt - v 2 )\ = \n ■ ((xi - x 2 ) - t*(vi - v 2 )) | < 3e ■ 
This tells us exactly that vi — v 2 belongs to the cylinder of axis X\ — x 2 and radius 3s /S. 

The lemma is proved. □ 

The second geometrical lemma requires the use of notation coming from scattering theory, introduced 
in Chapter 3: it states that if two points z\, z 2 in R 2d are in a post-collisional configuration and if v\ 
or v 2 belong to a cylinder as in Lemma 9.1.3, then the pre-image z 2 of z 2 through the scattering 
operator belongs to a small set of R 2d . 

Lemma 9.1.4- — Consider two parameters p <C R and i)Cl, and (y,w) € R d x Br. For any v\ 

in Br, define 

AA*(y,w,p)(ui) := {{v,v 2 ) G S? -1 x B R /(v 2 - Vl ) ■ v > n, 

vl e w + K(y,p) orv 2 e w + K(y,p)} , 
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where (v*,v*,V2) = cr Q 1 ""l, ^2) with the notations of Chapter 3. Then 

\M*(y,w,p)(v 1 )\<C($,R,r ] ) P d - 1 
where the constant depends on <& through the L°° norm of the cross-section b defined in Chapter 3. 

Proof. — Denote by r = \vi — U2I = \v* — v^l, and by u> the deflection angle. The formula (3.2.2) 
shows that, as u> varies in Sj -1 , v\ and v'2 range over a sphere of diameter r. 

The solid angle of the intersection of such a sphere with the cylinder w + K(y, p) is less than 

Gamin (l, (£) 

which implies that 

\{(w,v 2 )/vl ew + K{y,p) or v* 2 Ew + K{y,p)}\ < C d J r d ~ x min ^1, ^ dr 

< C d Rp d - x 

According to Chapter 3, the change of variables (u, V1—V2) i-> (u>, v\ — V2) is a Lipschitz diffeomorphism 
away from v ■ (v-y — v 2 ) — 0. We therefore get the expected estimate. □ 

Remark 9.1.5. — Note that those two lemmas consist in eliminating sets in the velocity variables 
and deflection angles only, and do not concern the position variables. 

9.1.3. Proof of Proposition 9.1.1. — We fix a good configuration Z k £ Qk(£o), an d we consider 
a configuration Z k G R 2dfc , with the same velocities as Z k , and neighboring positions: \X k — X k \ < a. 
In particular we notice that for all r > and all i ^ j, 

(9.1.7) \xi — Xj — r(v~i — Vj)\ > \xi — Xj — r(u, — Vj)\ — 2a > e /2 

since a -C £o- This implies that Z k £ Q k (eo/2). Next we consider an additional particle (x k +ev, v k +i) 
and we shall separate the analysis into two parts, depending on whether the situation is pre-collisional 
(meaning v ■ (v k +i — v k ) < 0) or post-collisional (meaning v ■ (v k +i — v k ) > 0). 

9.1.3.1. The pre-collisional case. — We assume that 

v ■ (v k +i - v k ) < , 

meaning that (x k + eu, v) and z k form a pre-collisional pair. In particular we have for all times r > 
and all e > 

j (x k +su- v k+1 r) - (x k - v k T) \>s. 
Furthermore up to excluding the ball B v (v k ) in the set of admissible v k +\, we may assume that 

\vk+\ -Vk\>V- 

Under that assumption we have for all r > 5 and all e > sufficiently small, 

j (xk +ev- v k+1 r) - (x k - v k r) \ > r\v k+1 -v k \-e 

>Sn-e> £ /2 • 

Furthermore we know that Z k belongs to Gk(£o/2) thanks to (9.1.7). 

Now let j e [1, k — 1] be given. According to Lemma 9.1.3, we find that for any v k+ i belonging to the 
set Br \ K(xj — x k ,QRa/eo + 6Rso/S), we have 

Vr > , I (x k + ev - v k+1 r) - (xj - vjt) \ > e, 
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and 

Vr > S , | (x k + ev - v k+1 r) - (x } - vjt) \ > e ■ 

Notice that 

\B R n K(xj - x k ,6Ra/e + 6Re /8) < C(R) ( — + ^ ) 
I \e o J 

Denning Mr(Z k ) := K(xj — x k , 6Ra/eo + 6Rso/6) and 
j<k-i 

B k (Z k ) -St 1 x (B n (v k )UM-(Z k )) 

we find that 

B- k (Z k )\<C{R)^ + (^) d " + (^) d ") 
and (9.1.3) and (9.1.4) hold as soon as (v, v k+ i) £ B k (Z k ). 

9.1.3.2. The post-collisional case. — We now assume that 

v ■ {v k +i -v k )>0. 

Next let us define 

(9.1.8) C(Z k ) := { (v, v k+1 ) eSixB R ,v (v k+1 - v k ) < rj} , 
which satisfies 

\C(Z k )\ < C(R)r 1 d - 1 . 

Choosing (v, v k +i) € (Sf x Br) \ C(Z k ) ensures that the cross-section is well defined (see Defini- 
tion 3.3.3), and that the scattering time t e is of order C(R,rj)s by Proposition 3.2.1. 

Considering the formulas (3.2.2) expressing (zf.*, zf+i) in terms of (z k , (x k + ev, v k+ {j), we know that 

(9.1.9) \x k -x%*\<e + Rt e <C(R,r])e and \(x k + ev) - x%* +1 \ < e + Rt e < C(R,t])e . 

Note that due to (9.1.7), all particles Xj with j < k — 1 arc at a distance at least eo/2 — e 3> e of the 
particles x k and x k +ev. Since they have bounded velocities, they cannot enter the protection spheres 
of these post-collisional particles during the interaction time t e , provided that e is small enough: 

Rt £ < e . 

Since the dynamics of the particles j < k— 1 is not affected by the scattering, we get that Z k *_ 1 belongs 
to £ fe _i(£ /2): 

(9.1.10) Vr > , V(i, j) e[l,k- l] 2 with ijtj, \xf - xf - t«* - vf)\ > e /2 . 

The pair (zf* , z k * +l ) is a prc-collisional pair by definition, so we know that for all r > 0, 

1(4* - tv%*) - (x%* +1 - rvf* +1 )\ >e. 
Moreover we have — t> £ *| = \v k — v k+ i\ > rj, so as in the pre-collisional case above we have 

Vr > 5 , \x%* - x s * - t(vI* - v s k * +1 )\ > V 6-e>e , 
for e sufficiently small, since e <C rjd. 
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Next for j < k — 1 we have for e sufficiently small, recalling that the uniform, rectilinear motion of the 
center of mass as described in (3.1.3), 

\xj* — Xj\ < \xj* — Xj | + \xj — xj\ < Rt e + a < 2a 

14* - x k \ < \\x%* -x% +1 \ + \\{xf + 4*+i) - (x k +x k+ i)\ + ^\(x k +x k+1 ) - 2x k \ 

< Rt e + e + a < 2a 

~ x k\ < ^14* -4*+i I + ^1(4* + 4*+i ) - (x k + x k+1 )\ + ^\(x k + x k+1 ) - 2x k \ 

< Rt £ + e + a < 2a . 

By Lemma 9.1.3, provided and v^* +1 do not belong to 

vj + K(xj - x k , 12Ra/e Q + Y2Rs /5) n B R , 



we get since vf* = v 



as well as 



Vr>0, \ x %* - x c* -t{v%* -v?)\>s, 
and \ x z* - x » - T {vi* +1 -v?)\>e 

Vr > 5/2 , 14* - 4* - 44* - «f )l > £o/2 , 
and - - r(C+i - wf )| ^ £ °/ 2 ■ 

A/"*(xj - 5 fe ,%,p)(ui) := {(v,v 2 ) S S^ 1 x B R / {v 2 - Vi) ■ v > jj , 

«i e w + if (y, p) or i>2 € w + p)} 

Lemma 9.1.4 bounds from the above the size of the set J\f*((xj — x k , Vj,p) of all (v, v k +i) € (Sf x Br) \ 
C(Z k ) such that vf.* or vf.* +1 belongs to Vj + K(xj — x k , p)- We let p = 12i?a/e + 12Rso/S, and define 

M + (Z k ) := |J JV*{xj -x k ,Vj,p) 
j<k-i 

and 

£+(Z fe ) := C(Z fc ) U (S^ 1 x X+(Z fe )) , 
where the set C(Z k ) is introduced in (9.1.8). By Lemma 9.1.4, 

B+(Z k ) < C{R) V d -' + C(<P, V ,R)(- + ^y- 1 

So 

and (9.1.5) and (9.1.6) hold as soon as (u, v) £ B^(Z k ). The proposition is proved. □ 

Note that, in order to prove that pathological sets have vanishing measure as e — > 0, we have to 
choose T) small enough, and then a and £o even smaller in order that (9.1.1) is satisfied and that (9.1.2) 
is small. Note that, if we want to get a rate of convergence, we need to have more precise bounds on 
the cross-section b in terms of the truncation parameters R and rj. 



9.2. Truncation of the Boltzmann and BBGKY pseudo-trajectories 

In this paragraph we show that the previous geometrical result may be used to define approximate 
Boltzmann and BBGKY observables, corresponding to non pathological pseudo-trajectories. We then 
expect to be able to compare these approximate observables, which will be done in the next chapter. 
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9.2.1. Initialization. — The first step consists in preparing the initial configuration Z s so that it 
is a good configuration. We define 

A s (e ) := \Z S £ R ds x B S R / inf \x e - xj\ > s \ , 

l l<£<j<s J 

and we shall assume from now on that Z s belongs to A s (e )- We also define for convenience 

Af ( £o ) := {X s £ K ds / inf \x t - Xj \ > e Q \ . 

Proposition 9.2.1. — For all X s £ A^-(e ), there is a subset M S (X S ) o/R ds such that 

\M s{ X s) \ KCfrs)^ 1 + 

and defining V s := \z s £ A s (e ) / V s £ M S (X S )^, then 

Vr>0, 1 Pi oH s (t)e1 Ps oS s (t) 
(9 ' 2 ' H) Vr><5, l P ,oS s (r)El Ps oS s (r)ol eiM . 

denoting abusively by TLa the operator of multiplication by the indicator of A. 

Proof. — The proof is very similar to the arguments of the previous paragraph. For any Z s in A s (e ), 
we apply Lemma 9.1.3 which shows that outside a small measure set M S {X S ) C R ds of veloci- 
ties {v\, . . . , v s ), with 

\M s (X s )\<C(R)C 2 s ^ + £ f^ 1 , 

the backward nonlinear flow is actually the free flow and the particles remain at a distance larger 
than e to one another for all times: 

Vt>0, Wj=£' £{1,...,s}, \(x t - vtr) - (x t , - v t ,r)\ >e, 

and that 

Vr><5, W^£' £{l,...,s}, \{x e -v e T)-(x e >-v ei T)\>e . 

By construction, M S (X S ) depends continuously on X s . Therefore, defining V s := {Z s £ A s (e ) / V s £ 
M S (X S )} gives the result. □ 

9.2.2. Approximation of the Boltzmann functional. — We recall that we consider a family of 
initial data F a — (f^) satisfying 

|||-Fb|||o^o,w = snp sup (exp(/3 E(Z s ) + fi s)f^ s) {Z s )) < +oo 

sGN Z s 

and after the reductions of Chapter 8, the observable we are interested in is the following: 

i° s £- 5 (t, J,M)(x s ) := / vs(v s ) ( s a (t-t 1 )<^r i s a+ i(ti -*2)C a °^; a+ , 

(9.2.12) J JT n At) 



ni2 



■ ■ ■ S s +„(t„ — t n +i)1L\Eo(z s+n )\<R 2 fo S+U '' ] dT n dV s , 
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By Proposition 9.2.1, up to an error term of order C(R)s 2 (e/en + So/S) d 1 , we can assume that the 
initial configuration Z s is a good configuration, meaning that 



j, M)(x s ) = f Vs (v s ) I s s (t ti)c a ° i ^r i s a+ i(ti - t 2 )c°4r™ 2 

JB R \Ms(X s ) JT n ,s(t) 

■ ■ ■^°s+^'-l,s+n^s+n{tn - tn+l)^\E (Z s+rl )\<R 2 fa^^ dT n dV s 



+0(C(R) 



where 



(C a s \+T.f {s+1) )(Z s ) = f ((v s+1 - v m ) ■ v s+1 )-.f {s+1 \Z s ,x m ,v s+1 )dv s+l dv s+1 and 

(Cs't+if {s+1) )( z s) = / {(vs+i -v m ) ■ iy s +i)+f {s+1) (z-i, . . . ,x m ,v^, . . . , z s ,x m ,v* +1 ) du s+1 dv s+1 . 

Jsf _1 xR d 

Now let us introduce some notation which we shall be using constantly from now on: given Z s e A s (e ), 
we call Zg(r) the position of the backward free flow initiated from Z s , at time t\ < r < t. Then 
given ji e {+,—}, mi G a deflection angle w s+ i and a velocity v s+ i we call Z° +1 (t) the 

position at time t 2 < t < t\ of the Boltzmann pseudo-trajectory initiated by the adjunction of the 
particle (v s +i, v s+i) to the particle (<i) (which is simply free-flow in the pre-collisional case ji = — , 
and free-flow after scattering of particles z° mi {ti) and (v s +i,v s+ i) in the post-collisional case ji = +). 

Similarly by induction given Z s £ A s (e ), T,J and M we denote for each 1 < k < n by Z® +k (r) 
the position at time tk+i < t < t k of the pseudo-trajectory initiated by the adjunction of the parti- 
cle {v s +k, v s+ u) to the particle z^ rik (tk) (which is simply free-flow in the pre-collisional case jk = — , 
and free-flow after scattering of particles z^ k (tk) and (v s+ k,v s+ k) in the post-collisional case jk = +)■ 

Notice that r i->- Z° +/c (t) is pointwise right-continuous on [0, tk]- 

With this notation, the functional Ig'n' S may be reformulated as 

I°;*> s (t,J,M)(X s )= f dV s ip s {V s ) ( dT n f dv s+1 dv s+1 ((v s+1 - v^ih) ■ v s+1 )+ 

JB R \M a (X s ) JT n , s {t) JSf^xBn 

... j dlS s+n dv s+n ((v s+n - V° mn {t n ) ■ V s+n ) + \ Eo{z o ( ))<R2fo +n \ Z s+n(°)) 

JSf- 1 xB a 



Let a, e , V -C 1 be such that 

a <C e <C 775 . 

According to Proposition 9.1.1, for any good configuration Z s+ k-i € ^d(s+k-i) j we can define a set 

c S s+ fe_i(Z s+ fe_i) := (S^ 1 x B R ) \ B s+ k-\{Z s+ k-i) , 

such that good configurations Z s+k -i = (X s+ k-i, V s+ k-i) with \X s+k _i — X s+k -i\ < Ca are stable 
by adjunction of a collisional particle {vk+s, Vk+s) <= c ^i+fe-i^s+fe-i)- 

We further notice that thanks to Remark 9.1.2, if the adjoined pair (v s+ k 1 v s+ k) belongs to the 
set c B s+ k-i{Z f} s+k _ 1 (t k )) with Z° s+k _ x (t k ) € Qs+k-\{^), then Z° s+k (t k +i) belongs to G s+ k(e )- 
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(9.2.13) 



As a consequence we may define recursively the approximate Boltzmann functional 
J°£' 5 (t,J,M)(X s )= f dV s <p s (V s ) f dT n 

JB R \Ms(X s ) JTnAt) 

/ dv s+1 dv s+1 (v s+1 - (ii) • v s+1 ) n 

... j du s+n dv s+n (v s+n - v° mn (t n ) ■ v s +n)j n 

J=B. +n _ 1 »(Zj +n _ 1 (t n )) 

x ^E„(Z° s+ JO))<R 2 fo +U) ( Z s+n( Q )) ■ 

The following result is an immediate consequence of Proposition 9.1.1 
Proposition 9.2.2. — Let a,eo,n <C 1 satisfying (9.1.1). Then, 
1a s(£0 )(/ s °^ 5 ~ J°;«' S )(t, J,M)\ < (s + nf (o(^ + C W M(J + f)^) ll^o||o,ft, w 



9.2.3. Approximation of the BBGKY functional. — We recall that after the reductions of 
Chapter 8, the functional we are interested in is 



i%(t,j,M)(x a ) := / Mv s ) I u s (t-h)ci^u s+1 ( tl -t 2 )C£?; s+2 

J JT n ,s(t) 

■ ■ ■^i+n-l,s+n^s+n(tn ~ tn+l)^\E e (Z s+n )\<R 2 fNfi n ^dT n dV s , 

where F Nfi = (7w,o)i<«<iv satisfies 

|||-FW,o||U,/3 ,m = supsup(exp(/3 S e (Z s ) +^ s)/^ s) (^ s )) < +°o • 
Thanks to Proposition 9.2.1, we have 

I^ s (t,J,M)(X s )= f Vs{V s ) I S a (t-t 1 )l j ,. (eo) C^ 1 1 H 8+1 (t 1 -t 2 )C^ +2 

Jb r \M 3 (X s ) JT n ,s(t) 

■ • •Cf+^!_ n ljS+ „H s+ „(t n - t n+ i)l Ee ( Zs+n (o))<R 2 fNfi n ^dT n dV s 



Then using the notation introduced in the previous paragraph for the Boltzmann pseudo-trajectory, 
let us define the approximate functionals 



J™(t,J,M)(X a ):= [ cp s (V s ) [ S a (t-ti)l j ,. (eo) Cf^ 1 1 H a+ i(t 1 -t 2 ) 

JBr\Ms{X s ) JT n ,5(t) 

■ ■ -Ci+n-l,s+nH-s+n(tn - t n +l) ^-E e (Z s+n (0))<R 2 fo*^ dT n dV s , 
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where 

(C;+T-i s+k 9 (s+k) ){Z s+k-i) := (N - a - k + l)^" 1 / {v s+k ■ (v s+k - v mh (t k )))- 

xg {s+h) {-,x mk (t k ) + eu s+k ,v s+k (t k )) Yl h(xj-xm k )(t k )-ev s+k \>edl> s+k dv s+k 

l<j<s+k-l 

(C+ + r-i. s+fe .9 (s+fe) )(^+fe-i) :=(N-s-k+ l)^" 1 / {v s+k ■ (v s+k - v mk {t k )))+ 

xg( s+k \. ..,x* mk ,v* mk ,..., x* +k ,v* s+k ) Y[ %\(x J -x mk )(t k )- ev .+ k \>e dv s+k dv s+k . 

l<j<s+k-l 

denoting as previously by {x* mk , v* mk , x* s+k ,v* s+k ) the pre-image of (x mk , v mk (t k ),x mk {t k )+su s+k , v s+k (t k )) 
by the scattering a e . 

As in the Boltzmann case described above, the following result is an immediate consequence of Propo- 
sition 9.1.1. 

Proposition 9.2.3. — Let a,So,r] <C 1 satisfying (9.1.1). Then, for e sufficiently small, 
l Aa{so) {l^-J^)(t,J,M)\ < (s + nf ^( J RV- 1 +C(i?,$, ?? )(^ + |)^^ \\\F Nfi \\\ E ^ . 

The functional can be written in terms of pseudo-trajectories, as in (9.2.13). Let us therefore 
introduce some notation which we shall be using constantly from now on: given Z s G A s (eo), we 
call Z®(t) the position of the backward free flow initiated from Z s , at time t\ < r < t. Then 
given ji e {+, — }, mi € [1, s], an angle v s+ \ (or cquivalently a position x s+ \ = x m ^ (t\) + ev s+ i) and 
a velocity v s+ i we call Zf +1 (r) the position at time t<i < t < t\ of the BBGKY pseudo-trajectory 
initiated by the adjunction of the particle z s+i to the particle z^iti)- 

Similarly by induction given Z s € A s (sq), T, J and M we denote for each 1 < k < n by Z| +fe (r) the 
position at time t fe+1 < r < t k of the BBGKY pseudo-trajectory initiated by the adjunction of the 
particle z s+k to the particle z mk (t k ). We have 

j£?(t, J, M)(X S ) = ( ^ N ~ S) l^ n(d - 1] ! dV s ip s (V s ) I dT n 

(W-s-n)! Jb r \m s {x s ) JT n , s {t) 

/ du s+1 dv s+1 (f s+1 ■ (V s+1 - !) mi (*l)))ii T] ^\(x j -x m .)(t 1 )-e Vs+1 \>e 

(9.2.14) Mmi 

(v 

s-\-n 

X II l|(^-^„)(t„)-£^ + „|>£lBa2 3+ „(0))<fl2/^o" ) (^+n(0)) • 

l<j<s+n-l 

Thanks to Propositions 9.2.2 and 9.2.3 the proof of Theorem 4 reduces to the proof of the convergence 
to zero of — Jg'n' S - This is the object of the next chapter. 



CHAPTER 10 



CONVERGENCE PROOF 



In this chapter we conclude the proof of Theorem 4 by proving that J^ 5 — J^n' S § oes to zero in the 
Boltzmann-Grad limit, with the notation of the previous chapter, namely (9.2.13) and (9.2.14). The 
main difficulty lies in the fact that in contrast to the Boltzmann situation, collisions in the BBGKY 
configuration are not pointwise in space, nor in time. At each collision time tk a small error is therefore 
introduced, which needs to be controled. 

We recall that, as in the previous chapter, we consider dynamics 

- involving only a finite number s + nof particles, 

- with bounded energies (at most R 2 » 1), 

- such that the n additional particles are adjoined through binary collisions at times separated at least 
by 5 < 1. 

The additional truncation parameters a,£o,?7 <C 1 satisfy (9.1.1). 
10.1. Proximity of Boltzmann and BBGKY trajectories 

This paragraph is devoted to the proof, by induction, that the BBGKY and Boltzmann pseudo- 
trajectories remain close for all times, in particular that there is no recollision for the BBGKY dynamics. 

We recall that the notation Z®(t) and Z k (t) were defined in Paragraphs 9.2.2 and 9.2.3 respectively. 

Lemma 10.1.1. — Fix T € T n ,s{t), J , and M and given Z s in A s (sq), consider for all i G {1, . . . n}, 
an impact parameter v s+ i and a velocity v s+ i such that {v s+ i,v s+ i) ^ S s +j_i(i?g +i _ 1 (tj)). Then, for e 
sufficiently small, for all i e [l,n], and all k < s + i, 

(10.1.1) \x%(ti+i) - x° k (t i+ i)\ < Cei and v k {t l+l ) = v° k (U +1 ) , 
where the constant C depends only on r],R,5. 

Proof. — We proceed by induction on i, the index of the time variables t i+ i for < i < n. 

We first notice that by construction, Z s (ti) — Z®(ti) = 0, so (10.1.1) holds for i = 0. The initial 
configuration being a good configuration, we indeed know that there is no possible recollision. 

Now let i e [1, n] be fixed, and assume that for all I < i 

(10.1.2) Vfc <* + *-!, \x%(t t ) - x° k (t t )\ <Cs(£-l) and v k (t e ) = v° k (t e ) . 
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Let us prove that (10.1.2) holds for I = i + 1. We shall consider two cases depending on whether the 
particle adjoined at time U is pre-collisional or post-collisional. 



• As usual, the case of pre-collisional velocities (v s +i,v mi (U)) at time U is the most simple to handle. 
We indeed have Vr G [ti+i,t»] 

\Jk<s + i, x&t) = x° k (U) + (t - U)vUU) , v k (r) = v° k (ti) , 

x "+ t ( T ) = x° m .{ti) + (r - U)v s+t , v s+t (r) = v s+l . 

Now let us study the BBGKY trajectory. We recall that the particle is adjoined in such a way 
that {v s+i , v s+i ) belongs to c B s+i _i{Z® +i _ 1 (t i )). Provided that e is sufficiently small, by the induction 
assumption (10.1.2), we have 

\X e a+i _AU) - X a ° +i _i(*i)l < Ce(i - 1) < a. 

Since Z^ +i _ 1 (t i ) belongs to Gs+i-i(eo) (see Paragraph 9.2.2), we can apply Proposition 9.1.1 which 
implies that backwards in time, there is free flow for Z £ s+i . In particular, 

Vfc<s + z, x k {r) = x k (U) + (t - ti)v k (U) , v k (r) = v k (U) , 

x s +i(r) = x m .(U) + sv s+ i + (t - U)v s+i , v s+ i(t) = v s+i . 

We therefore obtain 

(10.1.3) Vk<s + i, W€[t i+1 ,ti], v k (r) - v° k (r) = v k {ti) - v° k (ti) = , 
and 

(10.1.4) Vfc<s + «, Vr e [U+i,U]> \x k {r) - x k {r)\ < Cs{i - 1) + e . 



Vij+i <t<U 



• The case of post-collisional velocities is a little more complicated since there is a (small) time interval 
during which interaction occurs. 

Let us start by describing the Boltzmann flow. By definition of the post-collisional configuration, we 
know that the following identifies hold: 

' (i,."!+.)( T ) = with W+i.«m 4 (*i).^+i) : = <T0 1 K-H. u m 1 (*i)> w «+i) 

<M - a&,(*i) + (r - , x° +4 (r) = x° +4 (t 4 ) + (r - U)v* s+i 

Vj £ {m„ S + 1} , «?(t) = «?(ti) , x°(r) = arPfe) + (r - , 
where exo denotes the scattering operator defined in Definition 3.2.1 in Chapter 3. 

First, by Proposition 9.1.1, we know that for j ^ {mi, s + i} and W e [t i+ i 7 ti], 

Xj (r) = a:j (tj) + (r - t^u,- , Wj (t) = Vj (U) , 
so that by the induction assumption (10.1.2) we obtain 

Vjt{mi,8 + i},VT€[ti+i,U], \x j (T)-x° j (T)\ = \x j (t i )-x° j (t i )\<Cs 

(10 - L5) J f \ Of ^ 

and Vj(T) = Vj (r) . 



We now have to focus on the pair (s + i,rm). According to Chapter 3, the relative velocity evolves 
under the nonlinear dynamics on a time interval [U — t £ ,tj\ with t e < C(ij,R)e (recalling that by 
construction, the relative velocity \v s+ i — v m .(ti)\ is bounded from above by R and from below by 77, 
and that the impact parameter is also bounded from below by 77). Then, for all r G [i»+i,tt — t e ], 

(10.1.6) v s+l (r) = v* s+l = v° s+i (r) , v mi {r) = v* mi (U) = <(i.) = v° mi (r) . 
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In particular. 



(10.1.7) 




S-\-l 



(t i+1 ) and v mi (t i+ i) 



« mi (*i+i)- 



The conservation of total momentum as in Paragraph 9.1.3.2 shows that 



= a^fe) - + I < - 1) + | • 




(10.1.8) \x e m% {U - t e ) - x mi (U - t e )\ < Cei and \x e s+l (U - t e ) - x° s+l (U - t e )\ < Cei 



provided that C is chosen sufficiently large (depending on R and rf). 

Now let us apply Proposition 9.1.1, which implies that for all r G — t e ] the backward in time 

evolution of the two particles x e s+i (ti — t e ) and x e m .{ti — t e ), is that of free flow: we have therefore, 
using (10.1.6), 

x s+i(^+l) — ^s+i^i+l) = x s+i(ti ~ ^e) — x s+i(^i — ^e) ■ 

From (10.1.8) we therefore deduce that the induction assumption is satisfied at time step ij+i, and the 
proposition is proved. □ 

Note that, by construction, 



Note also that the indicator functions are identically equal to 1 for good configurations. We therefore 
have the following 




Corollary 10.1.2. — Under the assumptions of Lemma 10.1.1, the functional J^(t,J,M) defined 
in (9.2.14) may be written as follows: 



■ffift, J, M){X S ) = (JV S)! g***- 1 ) / dV s ip s (V s ) f dT n 

{N-s-n)\ Jb r \m s (x s ) JTnAt) 




x ^E e (Z 3+n (0))<R 2 ^ 



z s+n (o)eg s+n (e /2)/#,r ; (^ + „(o)). 
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10.2. End of the proof of Theorem 4 

The end of the proof of Theorem 4 consists in estimating the error terms in — J a 'n' S coming 
essentially from the micro-translations described in the previous paragraph and from the initial data. 

10.2.1. Error coming from the initial data. — Let us replace the initial data in by that 
of the Boltzmann hierarchy, defining: 

J* s (t,J,M)(X a )= ( N ~ s)l e n «-» [ dV aVa {V a ) [ dT n 

(jy-s-n)\ Jb r \m s (x s ) JT n ,s{t) 



I dv a+1 dv a+1 (v a+1 ■ (v s+1 - v mi (t 1 ))) jl 

J°B.{.Z°{ti)) 



■I 

^«=B s+n _ 1 (Z° +n _ 1 ( 



dv a+n dv a+n {v s+n ■ (v a+n - v mn (t n )))j n 
x ^E (z s+n (o))<R 2 ^zi + jo)eg s+ri (s /2)fo S+n) ( z s+n(0)) ■ 



Lemma 10.2.1. — In the Boltzmann- Grad scaling Ne d 1 = 1. 

I^Af (e )(J^n -^n)(*i J, M )( X s)\ < C — ||^ s ||i,oo( R i«) || 1a s+ „( £ „/2) (/jv,o" ~ fo" " ) || L oo (R 2d(.+n)) 

and in particular 

\TL A ?(e o ){T S ,n-Js,n)(t,J,M){X a )\^0 US E -> , 

uniformly int € [0, T] and X s e R ds . 

Proof. — We recall that by assumption, 1a s+ „(£ /2)(/v o™' — /o S+ ™' > ) § ocs to zero uniformy in Z s+n . 

By definition of the good sets Qk{c), the positions in the argument of Z^o™^ — fo S+ ^ satisfy the 
separation condition \xi — Xj\ > Eq/2 for i ^ j : 

JL e s+ „(£o/2)l/v,o -Jo ) - ll e s+ „( £ ,/2) J1 A J+ „( e „/2)H J v ! o -Jo )■ 
Furthermore, due to that separation condition, 

E £ (Z s+n ) = E (Z s+n ) = - \ v *\ 2 ■ 



2 

i=i 



So we can write 



(J^ 6 (t,J,M) - J^(t,J,M))(X s ) = ( N S) ' £ " (d - 1} / dV 6 <p.<y.) I dT n 

(N-s-n)\ JB R \M a (x s ) JT„At) 

/ dv a+1 dv s+1 (v s+1 ■ (v a+ i - v mi (t 1 ))) jl 

J°B a (Z°( tl )) 

... j dv a+n dv a+n {v a+n ■ (v s+n - v mn (t n ))) jn 

^<=B s+n _ 1 (Z° +n _ 1 (t n )) 



v | -n ( f(s+n) As+n)-, 

X lE e (Zf +n (0))<iJ2lA s+n ( £o /2)(/jv,0 -JO )> 



and we find directly that 



^ C ~j H^s||L oo (R ds )||lA s+T1 (£ /2)(/iV,0 n ~ /<)* " )|| L oo( R 2d( s + n 
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The result is proved. □ 

10.2.2. Error coming from the prefactors in the collision operators. — As e — > in the 

Boltzmann-Grad scaling, we have 



(N-s)\ ^ 



£ n(d-l) _^ 1 _ 



(N-s-n)\ 
Defining 

J*£(t, J,M){X S ) = [ dV sVs (V s ) f dT n 

JB R \M s (X a ) JT n , s {t) 

I dv s+1 dv s+ i (v s+ i ■ (v s+1 - v mi (h))) jl 

(10.2.9) J°B.(zo( tl )) 



... dv s+n dv s+n {v s+n ■ (V s+n - V mn (t n )))j 

7c Bs+n _ l( zO +n _ 1 (t n )) 



x ^E (z s+n (Q))<K^z- +n (o)eg s+7l (e /2)fa S+n) (Z s+n (0)) , 
we have the following obvious convergence. 

Lemma 10.2.2. — In the Boltzmann-Grad scaling 7Ve d_1 = 1. 

|iA f (Ji^-jf;f)M,M)(^ 

10.2.3. Conclusion. — We can now compare the definition (9.2.13) of J®fi' 5 {t,J,M) 
J°£> s (t,J,M)(X s ) = [ dV s <p s (V s ) [ dT n 

JB R \M S (X S ) JT n , 8 {t) 

I dv s+l dv s+1 ({v s+1 - v^ih) ■ v s +i) n 

... I dv s+n dv s+n ((v s+n - v°(t n ) ■ is s+n ) jn 

y«B. +n _ii(Z» +B _ 1 (t„)) 
^E o (Z° s+rl (0))<R 2 fo +n \ Z s+n( Q )) ■ 

and the formulation (10.2.9) for the approximate BBGKY hierarchy 
J*'*(t,J,M)(X s )= [ dV s ip s (V s ) [ dT n 

JB R \Ms(X s ) JT n , s {t) 

I du s+1 dv s+ i (v a+1 ■ (v a+1 - v mi (ti))) n 

... I dv s+n dv s+n (v s+n ■ (v s+n - v mn (t n ))) jn 

^«B. +n _i(z; +n _ 1 (t„)) 
x ^E (z s+n m<R 2l zi + jo)eg s+n (e B /2).fo S+n \z s+ n{0)) . 

Lemma 10.1.1 implies that at time we have 

\X s+n (0) - X° s+n (0)\ < C(R, V )ns, and V s+n (0) = V s ° +n (0) . 

Provided that /q S+ ™' is continuous, we then obtain the expected convergence and this concludes the 
proof of Theorem 4. 
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Notice for instance that if /q is Lipschitz, then we have the following estimate. 
Proposition 10.2.1. — In the Boltzmann-Grad scaling Ne^ 1 = 1, 

|1a? (e„)(I° - J)(t,X s )\ < C(R, V )ne\\Vx s+n ti S+n) \\L~ , 
uniformly in t € [0, T] and X s G R ds . 



CHAPTER 11 



CONCLUDING REMARKS 



11.1. On the convergence rate : the particular case of hard spheres 

The method of proof described in this text gives actually a more precise statement regarding the con- 
vergence then Theorem 4 : gathering all the estimates together, we indeed obtain a rate of convergence. 
For general short-range potentials of interaction, this rate is not completely explicit since the constant 
arising in Proposition 9.1.1 in the estimate of pathological sets depends on $, rj and R through the 
cross-section b. 

This constant can be made explicit in particular cases, especially in the simple case of hard spheres 
since the deflection angle oj and the normal v coincide. 

Note that all the arguments can be reproduced in this case once the dynamics for fixed N and e is well- 
defined (without multiple collisions - see Chapter 2). Moreover there are important simplifications : 

- True marginals coincide with truncated marginals because of the non penetration condition. In 
particular, there is no more need of cluster expansions, which simplifies a little bit the existence 
proof, and Proposition 8.1.1 is no more relevant. 

- The scattering operator is completely explicit since lj = v in formulas (3.2.2). In particular the 
cross-section 



As mentioned above, this enables us to get a constant C(R, 77, 4>) = d in Lemma 9.1.4 that depends 
only on the dimension d. Since we use Lemma 9.1.4 with p = 12Ra/e + 12Re /S in the proof of 
Proposition 9.1.1, this gives an error 



Note that we can choose rj — C^f for C sufficiently large. 
- Collisions are pointwise and instantaneous t £ = 0, which makes the proof of Lemma 10.1.1 on the 
divergence of trajectories very easy. Indeed, the distance between the BBGKY and Boltzmann 
pseudo-trajectories increases at most of e at each collision! 

Let us then gather all the estimates together. We assume for the sake of simplicity that we start from 
an almost factorized initial data, i.e. a BBGKY initial data obtained from a tensor product by the 
conditioning process described in Chapter 7. 



b(vi - v 2 , w) = ((v! - v 2 ) ■ w)+ . 
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From the arguments of Proposition 8.2.1, 8.3.1 and 8.4.1, we find a first error term 



where Co depends only on the L°° norm of f - Then, from Proposition 9.1.1, we obtain as discussed 
above the error term 

e 2 <CoR d (- + £ °"" 



Finally, we have to take into account the error coming from the initial data, estimated in Chapter 7 
and Lemma 10.2.1: 

j^dnrpn 

e3 < C : — (s + n)e , 

n! 

the error coming from the prefactors of the collision operators 

e4 < Cqs{s + n)e 2 

and the error coming from the divergence of trajectories, which can be estimated if /o is Lipschitz as 
follows 

&b < C||/ ||wi,°o — — (s + n)e. 

Therefore, choosing 

n~Ci|log £ |, i? 2 ~C 2 |loge| 
for some sufficiently large constants C\ and C2, and 

S = £ (d-i)/(d+i)^ £Q ^ e d/(d+i) 

we find that the total error is smaller than Ce a for any a < (d — l)/(d + 1). 



11.2. On the time of validity of Theorems 2 and 4 

Let us first note that, for any fixed N, the BBGKY hierarchy has a global solution since it is equivalent 
to the Liouville equation in the phase space of dimension 2Nd, which is nothing else than a linear 
transport equation. The fact that we obtain a finite life span is therefore due to the functional 
spaces X e t p tft we consider. Belonging to such a functional space requires indeed a strong control on 
the high order correlations. The estimates we have written show actually (see Corollary 6.1.2) that 
the time of validity of Theorem 2 depends on /io , which measures the logarithmic growth of the initial 
marginals (that is the size of /o for factorized initial data) . 

An important point is that the time of convergence is exactly the time of existence. By definition of 
the functional spaces, we are indeed in a situation where the high order correlations can be neglected 
(see Proposition 7.2.1), so that we only have to study the dynamics of a finite system of particles. The 
term-by-term convergence relies then on geometrical properties of the transport in the whole space, 
which do not introduce any restriction on the time of convergence. 

A natural question is therefore to know whether or not it is possible to increase the time of existence 
and thus the time of convergence. The interpretation we have of the iterated Duhamel formula in terms 
of pseudo-trajectories, for which at each time step a new particle (independent from all the others) 
is added and interacts with one of the previous ones, shows that one cannot hope to prove by this 



11.3. MORE GENERAL POTENTIALS 
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method a convergence theorem for larger times than a fraction of the time at which each particle has 
undergone one collision. 

The first trivial remark is that after N — 1 iterations, independence does not hold any longer. Indeed, 
representing each collision by a strap between the two involved particles, we see that 

- either there is a chain connecting all the particles, 

or there are closed subchains, corresponding to pathological dynamics involving recollisions. 
It is actually known (see for instance [6]) that after cN time steps with 1/2 < c < 1, there is a change 
of phase and the appearance of a "giant component" in the set of N particles, meaning a set of a(c)N 
particles, with a(c) — >■ 1 as c — > 1 which have interacted either directly, or indirectly through other 
particles, so independence does not hold anymore, after N/2 time steps. 

In other words this means that, starting from such a system of N particles, we can expect to increase 
the time of convergence (and thus the time of existence) only if we can prove that, after a short time, 
particles go at infinity in different directions, and thus do not encounter each other any more, i.e. the 
dynamics reduces to free transport. 

If we want to establish the validity of the Boltzmann equation for longer times, we have therefore to 
start with systems of particles which contain much more particles, but with the same average density 
(in order that the collision cross-section remains bounded). The difficulties are then to prove that 

- the density of particles remains locally bounded, so that the asymptotics is still governed by the 
Boltzmann-Grad scaling; 

- the spatial dispersion creates some mixing mechanism, which implies that particles entering a colli- 
sion are always independent. 

Note that a simple way to get rid of the first problem is to consider periodic distributions of particles. 
In that case, the challenge is to understand how the dispersion associated to free transport could help 
for the propagation of chaos, which implies more or less to study the spatial decay of correlations. 

11.3. More general potentials 

A first natural extension to this work concerns the case of a compactly supported, repulsive potential, 
but no longer satisfying (1.2.1) of Assumption 1.2.1. As explained in Chapter 3, that assumption 
guarantees that the cross section is well defined everywhere, since the deflection angle is a one-to-one 
function of the impact parameter. If that is no longer satisfied, then one expects that additional 
decompositions are necessary, and resummation procedures need to be justified (see [36]). 

From a physical point of view it would be more interesting to study the case of long-range potentials. 
Then the cross section actually becomes singular, so a different notion of limit must be considered, 
possibly in the spirit of Alexandre and Villani [3]. One intermediate step, as in [15], would be to 
extend this work to the case when the support of the potential goes to infinity with the number of 
particles. Then one could try truncating the long-range potential and showing that the tail of the 
potential has very little effect in the convergence. 

Note that in the case when grazing collisions become predominant, then the Boltzmann equation 
should be replaced by the Landau equation, whose derivation is essentially open; a first result in that 
direction was obtained very recently by A. Bobylev, M. Pulvirenti and C. Saffirio in [4], where a time 
zero convergence result is established. 
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11.4. Other boundary conditions 

As it stands, our analysis is restricted to the whole space (namely Xn <G R dN ). It is indeed important 
that free flow corresponds to straight lines (see in particular Lemmas 9.1.3 and 9.1.4 as well more 
generally as the analysis of pathological trajectories in Chapter 9). 

It would be very interesting to generalize this work to more general geometries. A first step in that 
direction would be to study the case of periodic flows in Xn- The geometric lemmas must be adapted 
to that framework, and in particular it appears that a finite life span must a priori be given before the 
surgery of the collision integrals may be performed. The case of a general domain is again much more 
complicated, and results from the theory of billiards would probably need to be used. 
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B 8 R , ball of radius R centered at zero in R ds , 
page 17 

Br{x), ball of radius R centered at x in R d , 
page 32 

Br, ball of radius R centered at zero in R d , 
page 17 

Bk (Zk) a small set of angles and velocities of a par- 
ticle adjoined to Zk (or a neighboring configura- 
tion), leading to pathological trajectories, page 66 

b, cross-section, page 20 

C N , BBGKY hierarchy collision operator, page 29 

C°, Boltzmann hierarchy collision operator, 
page 29 

C s . s+m7 BBGKY collision operator involving m ad- 
ditional particles, page 28 

C° s+1 , Boltzmann collision operator, page 37 

T> N , artificial set on which the dynamics takes 
place, page 22 

D s (t), s-particle flow for hard spheres, page 10 

A m (X s ), m-particle cluster based on X s , page 26 

A s , set of well separated initial configurations, 
page 71 

, set of well separated initial positions, page 71 
da 1 ^ , surface measure on Y, s N (i,j), page 24 
da, surface measure on S s (xi), page 27 



dXij, d(j — i + l)-dimensional Lebesgue measure, 
page 22 

dZ iy j, 2d(j — i + l)-dimensional Lebesgue measure, 
page 22 

E(X s ,X n ), e-closure of X s in Xn, page 26 

E < i j 0> (X s ,X n ), e-closure of X s in Xn with a 
weak link at (io,jo), page 26 

E S (Z S ), s-particle Hamiltonian, page 34 

E (Z S ), s-particle free Hamiltonian, page 34 

(s) 

f N , marginal of order s of the iV-particle distri- 
bution function, page 21 

fffi, truncated marginal of order s of the TV- 
particle distribution function, page 22 

/( s ) , marginal of order s associated with the Boltz- 
mann hierarchy, page 39 

<fr e , rescaled potential, page 34 

Qk, set of good configurations of k particles, 
page 65 

H s (t), s-particle flow, page 28 

H(t), BBGKY hierarchy flow, page 29 

I v , observable (average with respect to momentum 
variables), page 52 

K{y,rj), cylinder of axis y and radius r\, page 67 
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Kd, volume of the unit ball in R d , page 31 

i/*'- 7 ', direction of Xi — Xj, page 24 

M S (X S ), good set of initial velocities associated 
with well separated positions, page 71 

p*, distance of minimal approach, page 15 

S s (t), s-particle free flow, page 29 

S^ -1 , unit sphere in R d , page 14 

S e (xi), sphere in R d of radius e, centered at x i} 
page 27 

cr e , scattering operator, page 16 

cr , Boltzmann scattering operator, page 16 

Sjv(i, j), boundary of the artificial set T> S N , page 24 

t e = nonlinear interaction time, page 15 

T n (t), set of collision times, page 61 

T n ,s(t), set of well-separated collision times, 
page 61 



-X"e, s ,/3 function space for the BBGKY marginals, 
page 34 

X £j/ 3 jAi function space for the BBGKY hierarchy, 
page 34 

X £!l g. M function space for the uniform existence to 
the BBGKY hierarchy, page 39 

Xq^ju function space for the Boltzmann hierarchy, 
page 34 

w, direction of the apse line, page 16 

r2jv, phase space for the Liouville equatoin, 
page 21 

Z^v, partition function, page 46 

| • | ejS ,^ norm for the BBGKY marginal of order s, 
page 33 

I ' \o,s./3 norm for the Boltzmann marginal of or- 
der s, page 33 

IMIe,^ norm for the BBGKY hierarchy, page 34 

II • II norm for the Boltzmann hierarchy, page 34 

II ' lk/3,M> n0rm in X £,/3,M' P a S e 39 



